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PREFACE 


A.  Target 

This  report  is  designed  for  YOU,  the  user,  and  is  so  written. 


B.  Notation  and  Usage 

The  numeral  'O'  (zero)  and  the  letter  'O'  (oh)  are  easily  confused. 
In  this  report,  the  character  'O'  is  always  the  letter  when  it  appears 
with  other  letters,  as  in  COM  or  ORD.  If  it  appears  with  numbers,  as  in 
3.04  or  1.0.  it  is  a  numeral.  In  a  few  places  where  misunderstanding 
may  occur,  the  symbol  '0'  is  used  for  the  letter  'oh'.  (This  is  standard 
BRL  usage.) 

In  this  report,  certain  names  appear  inside  single  quote  marks,  as 
'name2'.  When  this  happens,  you  are  to  understand  that  it  is  the  quan¬ 
tity  denoted  by  the  symbol  'name2'  that  is  being  referred  to,  not  the 
characters  themselves.  (This  is  contrary  to  standard  computer  usage, 
but  is  adopted  here  m  order  to  minimize  the  use  of  quote  marks.) 

In  standard  American  English  usage,  punctuation  marks  go  inside 
quote  marks.  In  this  report,  however,  the  punctuation  is  placed  outside 
the  quote  marks  to  decrease  the  likelihood  of  confusion. 


C .  Applicability  of  This  Report 

Since  late  1978,  BLAKE  has  printed  a  version  designator  on  every 
output  page.  This  report  is  specifically  intended  for  use  with  version 
20S.0.  Any  change  in  the  program,  no  matter  how  small,  will  be  signalled 
by  a  change  in  version  number. 

In  any  case,  the  nature  of  the  change  or  correction  will  be  nv  .ed 
in  comment  cards  placed  in  Subroutine  VRSION.  If  necessary,  chants  to 
this  manual  will  be  noted  there  also. 


II.  INTRODUCTION 

The  BLAKE  computer  program  is  a  general  equilibrium  thermodynamics 
program  derived  from  an  older  version  of  TIGER.  Although  it  is  applicable 
to  a  wide  range  of  chemical  equilibrium  calculations,  it  is  specifically 
intended  for  computing  the  properties  of  gun  propellants  at. chamber  con¬ 
ditions. 


7 


A.  Some  Background 


Among  the  first  chemical  equilibrium  programs  that  did  not  assume 
the  ideal  gas  equation  were  BKW  and  RUBY.  They  both  employed  the  Becker- 
Kistiakowsky-Wilson  equation  of  state  that  had  been  developed  earlier 
at  Los  Alamos  for  the  computation  of  p-V-T  properties  behind  detonation 
waves  in  condensed  explosives.  RUBY  apparently  was  difficult  to  use, 
Tequiring  considerable  user  supervision  and  intervention.  To  improve 
this  situation,  the  US  Army  Ordnance  Corps'  Ballistic  Research  Laborato¬ 
ries  (now  the  Ballistic  Research  Laboratory  of  the  US  Army  Armament 
Research  and  Development  Command)  contracted  with  the  Stanford  Research 
Institute  (now  Stanford  Research  International)  to  develop  an  entirely 
new  code  that  would  be  much  easier  to  use  and  that  would  permit  the 
ready  insertion  of  a  variety  of  non-ideal  equations  of  state  for  the 
gaseous  and  condensed  states.  The  TIGER  code  was  the  result.  Since 
then,  managerial  oversight  of  TIGER  has  resided  in  a  group  representing 
the  Army,  Navy,  and  Energy  Departments.  Under  its  guidance,  the  pro¬ 
gram  has  been  significantly  improved  and  its  capabilities  extended. 

Originally  TIGER  incorporated  only  two  gaseous  equations  of  state, 
the  ideal  and  the  BKW,  neither  of  which  is  applicable  to  gases  at  gun- 
chamber  conditions.  These  c-re  typically  pressures  up  to  700  MPa  and 
temperatures  between  1500  and  3000  K,  too  dense  for  the  ideal  gas  equa¬ 
tion  but  not  nearly  dense  enough  for  the  BKW.  A  more  appropriate  equa¬ 
tion  for  such  conditions  is  the  truncated  virial  equation, 

pV  =  nRT  [1  +  (n/V)B(T)  ♦  (n/V)2C(T)),  (1) 

whose  application  to  this  problem  was  first  described  by  Corner1. 

Corner  showed  how  the  constants  B(T)  and  C(Tj,  the  second  and  third 
virial  coefficients,  could  be  computed  once  the  intermolecular  potential 
was  known.  The  significance  of  this  result  is  that  it  bypasses  the  need 
for  measurements  on  the  chamber  gases,  a  complex,  mixture  of  some  10  to 
40  molecules,  free  atoms,  and  radicals,  at  elevated  temperatures  and 
pressures. 

A  code  embodying  Corner's  results  was  written  over  22  years  ago  in 
the  BRL  by  Leser^.  Results  obtained  with  it  were  summarized  in  a  com¬ 
prehensive  report  by  Baer  and  Bryson^  which  is  still  a  valuable  reference 
work. 

^J.  Comer ,  "Theory  of  the  Interior  Ballistics  of  Guns,"  J.  Wiley  &  Sons , 
NY  &  London  (1950). 

2 

T.  Leser s  "High  Speed  Digital  Calculations  of  the  Composition  and  Thermo¬ 
dynamic  Properties  of  Propellant  Gases  by  the  Brinkley  Method3"  BRL  R1023 
(August  1957)  (AD  156716) 

7 

P,  G.  Baer  and  K.  Bryson ,  "Tables  of  Computed  Thermodynamic  Properties 
of  Military  Gun  Propellants a"  BRL  MR  1358  (March  1961).  (AD  258644) 
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Leser's  program  was  written  in  assembly  language  for  ORDVAC.  Some 
years  later,  Baer  and  Wortman  produced  a  new  version  of  the  program  in 
a  high-level  language  called  FORAST.  The  Baer-Wortman  program  is  docu¬ 
mented  in  an  appendix  to  a  report  by  Stansbury  and  Shearer^*.  This  latter 
program  became  a  veritable  work -horse  within  the  BRL,  but  it  could  not 
be  exported  because  the  explosive  acceptance  of  FORTRAN  negated  the  pos¬ 
sibility  of  FORAST  compilers  being  written  elsewhere. 

When  the  use  of  FORAST  was  proscribed,  the  TIGER^  code  was  already 
available.  One  of  TIGER'S  many  outstanding  features  is  its  modular  design, 
which  enables  one  to  insert  other  gaseous  equations  of  state  by  changing 
only  one  subroutine.  Rather  than  translate  Baer-Wortman,  it  was  decided 
to  take  advantage  of  TIGER  and  thus  gain  the  benefits  of  TIGER'S  other 
attractive  features.  The  result  was  BLAKE. 

From  a  common  origin  in  1964  TIGER  and  BLAKE  have  proceeded  along 
separate  paths;  both  programs  have  changed  so  much  that  their  1979  ver¬ 
sions  cannot  use  input  cards  prepared  in  1964.  The  improvements  in 
TIGER  have  always  broadly  emphasized  the  computation  of  detonation  para¬ 
meters,  including  improved  convergence,  two-phase  systems,  and  partially- 
frozen  expansions.  Consequently  it  remains  the  code  of  choice  for  de¬ 
tonation  computations  for  condensed  explosives.  BLAKE  has  continued  to 
emphasize  gun  propellants;  in  its  present  form,  its  original  capabili¬ 
ties  for  computing  Chapman-Jouguet  parameters  and  Hugoniot  curves  have 
been  deleted. 


B ,  Comments  on  Gaseous  Equations  of  State 

The  truncated  virial  equation  is  nowadays  primarily  of  theoretical 
significance  only.  Practicing  chemical  engineers  have  at  their  disposal 
many  better  equations  (e.g.,  Benedict -Webb-Rubin,  Beattie-Bridgman) . 
These  equations  cannot  be  used  for  propellant-gas  computations  because 
there  is  no  feasible  way  of  finding  the  values  of  their  many  parameters; 
the  Benedict -Wcbb-Rubin  equation,  for  instance,  has  8-. 
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L.  Stansbury,  Jr.  and  R.  B .  Shearer,  " Tables  of  Computed  Thermo¬ 
dynamic  Properties  of  Selected  Liquid  Monopropellants,"  BRL  R1579  (April 
1972).  (AD  S21S26L) 

5 a)  W.  E.  Wiebenson,  Jr.,  W.  H.  Zwisler,  L.  B.  Seely ,  and  S.  R.  Brinkley,  Jr., 
" TIGER  Computer  Program  Documentation  (November  1968); 

b)  M.  Cowperthwaite  and  W.  H.  Zwisler,  "TIGER  Computer  Program  Documen¬ 
tation,"  SRI  Publication  No.  Z106  (January  1973). 
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The  significance  of  Corner's  work  is  that,  although  the  truncated 
virial  equation  is  far  from  the  best,  it  is  exceptionally  useful.  Corner 
used  the  Lennard-Jones  6,12  intermolecular  potential  function, 6 

4>(r)  =  4(e/kT)[(o/r)  -(o/r)  ]  (2) 

where 

<f>(r)  =  intermolecular  potential; 
e  =  depth  of  the  attractive  well; 
a  ~  distance  at  which  the  attraction  vanishes;  and 
r  =  distance  between  two  molecules, 
h  =  Boltzmann  constant,  and 
T  =  absolute  temperature 

The  exponent  6  is  dictated  by  theory;  the  exponent  12  was  chosen 
to  give  the  best  fit  for  second  virial  coefficients. 

A  drawback  of  the  Lennard-Jones  potential  is  that  it  is  inherently 
incapable  of  representing  the  intermolecular  potentials  of  either  non- 
spherical  or  polar  molecules,  especially  water.  The  Stockmayer  poten¬ 
tial7,  which  adds  the  term  -  y2G/r3  (y  =  dipole  moment,  and  G  =  a  function 
of  the  angles  that  specify  the  orientation  of  the  dipole  in  space)  to 
Eqn.  (1),  was  developed  to  overcome  this  objection. 

Unfortunately  the  situation  is  too  complicated  to  be  redeemed  by  a 
single  extra  term.  While  the  Stockmayer  potential  does  give  improved 
estimates  for  B(T),  it  does  not  improve  the  estimates  of  C(T). 

Numerically  the  contribution  of  the  third  coefficient  is  necessary  - 
for  loading  densities  greater  than  about  0.2  g/cm^,  but  theory  does  not 
appear  capable  of  giving  adequate  estimates  for  it.  For  this  reason, 

BLAKE  eschews  the  complications  inherent  in  making  extensive  theoretical 
estimates  of  C(T);  instead,  it  follows  an  earlier  re  „»>mmendation  of 
Hirschfelder's®  and  computes  the  third  virial  coefficient  with  the  as¬ 
sumption  that  the  molecules  are  hard  spheres  with  a  radius  equal  to 
0.81  times  the  Lennard-Jones  radius  o. 

For  a  thorough  discussion  of  gaseous  equations  of  state  and  their  re¬ 
lation  to  interwolecular  potentials 3  see  J.O.  Hir schf ‘elder 3  C.F.  Curtiss 3 
and  R.B.  Bird3  ", Molecular  Theory  of  Fluids 3"  John  Wiley  &  Sons 3  New  York 
(1954). 

? See  the  discussion  in  ref.  6. 

®See  ref.  63  p.  157  and  pp.  262-263. 
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This  decision  was  made  early  during  the  development  of  BLAKE.  In 
1977  a  coding  error  was  discovered  in  this  part  of  the  code;  in  retro¬ 
spect,  it  would  have  been  better  to  have  changed  the  model  from  hard 
spheres  to  square  wells  at  that  time,  but.... 

Corner's  approach  applied  strictly  only  to  A-A  and  B-B  interactions 
in  a  mixture  of  A  and  B;  A-B  interactions  were  ignored.  Although  it  may 
be  inconsistent  with  the  neglect  of  sophisticated  third-viria.1  computa¬ 
tions,  the  present  BLAKE  employs  mixing  rules  for  both  the  parameters 
of  the  intermolecular  potential  and  of  the  virial  coefficients.  They 
are: 


3  3  1/2  _3  1/2 

e. .  =  [a.  o.  1  '  o. .  [e . e  .1  , 

13  1  i  J  J  ij  i  3J 

(3) 

a..  =  0.5  (a.  +  a.); 
ij  i  3 

(4) 

B(T)  -  n.XiXjB..;  and 

(5) 

c(T)  =  yx.c _ 

(6) 

In  these  equations,  Xj  and  xj  are  the  mole  fractions  of  species  i  and  j 
(which  may  be  the  same) ,  is  the  virial  coefficient  for  the  mixture 
of  gaseous  species  i  and  j  (computed  using  the  parameters  j  and  j ; 
and  B(T)  and  C(T)  are  the  second  and  third  virial  coefficients  of  the 
total  mixture. 

Note  that  the  mixing  rule  for  the  third  virial  coefficient  is  simply 
a  linear  weighted  sum,  with  all  interactions  ignored.  This  is  consistent 
with  the  use  of  the  hard-sphere  model. 


Algebra  of  the  Truncated  Virial  Equation 


The  insertion  of  a  new  gaseous  equation  of  state  into  TIGER  requires 
the  evaluation  of  a  number  of  expressions  which  are  then  programmed  and 
put  into  Subroutine  STATEG. 


For  reference  purposes,  these  expressions  are  given  here  for  the 
truncated  virial  equation.  The  exposition  is  not  to  be  considered  a 
derivation. 


The  truncated  virial  equation  is 

p  =  (RT/V) [1  ♦  (B/V)  +  (C/V2)]  (7) 

for  1  mole  of  gas  of  volume  V.  Here  R  is  the  universal  gas  constant, 
and  H  und  C  are  the  second  and  third  virial  constants. 
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A  completely  general  equation  of  state  is 

P  =  nRT<t»p/Mo  ($) 

If  4>  =  1,  this  becomes  the  ideal  gas  equation  of  state;  in  general,  * 
is  a  measure  of  the  deviation  from  ideality. 

In  the  general  case,  B  and  C  are  functions  of  temperature,  but  here, 

C  will  be  assumed  independent  of  temperature. 

The  volume,  V,  is  related  to  the  reference  mass,  M  ,  and  the  density, 

P ,  by 


p  =  M  /V 
o 

The  following  quantities  will  also  be  needed: 


(9) 


X-  :  mole  fraction  of  i 
i 

n.  :  moles  of  i 

i 

n  :  total  moles  =  yn. 

-  i 

then. 


B  = 


r 


p  = 


IxjXjSij 

2 

y (n.n ./n  )  B.  . 

'-v  l  y  1  ij 

(10) 

y(n./n)3C. 

(11) 

nRTp<J>/M 

0 

(nRTp/Mo)[l  +  (Bnp/MQ)  +  (np/MQ)  C] 

(12) 

The  TIGER  formalism  requires  the  quantities  3P/3nj,  T^,  31n<*>/31nT, 
31n<5>/31np,  31n<t>/3nj,  ST^/SlnT,  3I\/31np,  3I\/3nj,  e,  and 

Some  of  these  have  already  been  defined,  the  others  are 


ri  =^P[(M0/RTp)C3p/3n.)-l]dp/p  (13) 

e  =/P(Mo/RTP)[p-CW3T)T3dp/p,  and  (14) 

e'T  =  f*  (M0T/Rp)(32p/3T2)dp/p  (15) 
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''".\c  origin  and  physical  meaning'/  of  all  of  these  quantities  are  in 
the  TIGP.R  documentation  and  will  not  be  duplicated  here. 

When  the  truncated  virial  equation  (with  constant  C)  is  substituted 
into  these  expressions,  tedious  aigebru  leads  to  the  following  equations: 


p  =  inRTp/M  }  +  (p/M  )  RT  Yn.n.B. . 

o  c  ^  n  11 

3  ..  3 

+  (p-M  )  KT  }n.  C. 
o  L  1  1 

2  2 

F.  =  (2p/M  )7n  B. .  +  (5/2) (p/M  )  n.  C. 
l  o  L  g  ij  w  M  o'  i  i 

2  3 

4  =  1  +  (p/M  n)Vn.n.B..  +  r.(p/M  n)  Yn.  C. 
cln$/3JnT  -  (p/d>)  ( (l/M^^n.n.-B. 

+  (2p/Mq  n)^n.  C.j 

Jlr.$/3n.  "  (p/M  n<5>)  2>n.B..  -  (l/n)YYn.n.B. . 

l  o  L  j  ij  ’U  i  1  ij 

♦  (P/M}(3n.2  -  £(n.  /n))Ci 
3Fi/31nT  *  (2p?/Mo)Yn.3B. ,/3T 
DTi/SlnP  »  Ti  +  (3/2)  (o/M  )2n.*C. 

3ri/3n.  «  (p/Mo)[2B.  .  ♦  (op/M^nA .0.  ] 

^•Cii  r‘*le  ^ronec{'ST  function,  defined  as 
ai-  5  0  if  i  /  j,  o..  =  l.j 

£  =  -lpT/Mo)nninj3B..m 

c’l  =  -(p/M0)nn.n,[2TdB  /dT 

2  2  J  2  ~ 

*  (T  /M  )d  B. ,/dT^J 
o  13  ■* 


(16) 

(17) 

(IS) 

(19) 

(20) 
(21) 
(22) 

(23) 

(24) 

(25) 

(26) 


I  want  to  thank  Arthur  Cohen  for  checking  the  algebra  (and  finding 
an  error),  and  Caledonia  lienry  for  double-checking  the  algebra  and  check¬ 
ing  the  programming  required  to  inject  these  equations  into  the  code. 
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m. 


COMPUTATIONAL  AND  P.XPHRSMHNTAL  ThSTS  OF 


'  BLAKF ' 


t 


- 
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A .  Computational  T e  st 

The  computational  testing  of  BLAKF  has  two  j..  poets.  The  first,  which 
is  readily  accomplished,  is  verifying  that  the  code  produces  the  correct 
results  for  thermodynamic  computations  on  ideal  gases.  One  could  do  this 
by  hand  calculations  but  the  ready  availability  of  several  wel 1-val i dated 
codes  makes  this  unnecessarily  laborious.  The  NASA -Leu is  code,  CEC'l®, 
was  taken  to  be  the  standard.  Computations  were  run  on  it  for  five  dif¬ 
ferent  propellant  formulations.  These  five  pro-pel  lanes  were  chosen  not 
as  representative,  current,  or  typical,  but  mainly  to  ensure  that  the 
temperature  range  from  2250  to  3800  K  was  covered,  and  that  at  least  one 
of  them  formed  a  condensed  phase.  Loading  densities  from  0.05  to  0.40  in 
intervals  of  0.05  were  used;  some  calculations  at  a  loading  density  of 
0.00  were  also  made.  All  calculations  used  the  same  enthalpies  of  forma¬ 
tion  as  input. 

A  code  developed  bv  the  Institut  fuer  Chem'e  dcr  Trcib-  und  Fxnlo- 
sivstoffe1  was  also  run. 

j } 

The  results  are  described  in  another  report  J ,  but  it  is  very  inac¬ 
cessible  so  they  will  be  repeated  here. 

Hie  five  compositions  are  shown  in  Table  i,  and  the  results  in  Table 
2.  The  concordance  is  excellent.  Additional  evidence  for  the  quality 
of  the  agreement  is  found  in  the  M15  calculations  where  8LAKH  and  Ci-.C'l 
both  predict  the  formation  of  liquid  A1  ->0-  with  mole  fractions  of  0.000124 
and  0.00012,  respectively. 

Hie  agreement  between  BLAKF  and  the  ICT  code  Is  shown  in  Table  3. 
ft  is  also  very  good,  but  not  quite  so  perfect  as  between  NASA- Lew  is  and 
BLAKF..  The  reason  likely  is  in  the  thermodynamic  data  selected  for  the 
products. 


"5.  Jordon  and  3.  J .  ”c3ride ,  "Computer  Pr-o.iv.ari  for  Computation  of 
Complex  Chemical  Equilibrium  Compositions _ NASA  CP-273  (1371). 

' J ? .  Volk  and  H.  Baikal t ,  "Application  of  the  Vivial  Equation  o  ' 
Ctate  in  Cala-'ating  Interior  Ballistic  Quantities."  Prove  ll-ar.  is  and 
Explosives  1.  7-PI  7 1976 ) 

't  i 

1 1 P.  Freedman,  "An  In  ter  comparison  of  West  Jovian  and  3S  Propellant 
Thermodynamic  Codes, "  in  Agenki  of  Presentations  at  the  25-23  April 
1973  Meeting  of  the  FRO  and  UP  Personnel  Associated  with  PEA  1060,  p. 
IB-V-119  (April  1973). 
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TABLE  2.  COMPARISON  OF  BLAKE  WITH  NASA/LEWIS 


Although  they  employ  different  fittings  for  them,  BLAKE  and  CEC71 
both  use  virtually  the  same  numbers  from  the  JANNAF  Tables*^  it  is  not 
clear  which  data  were  used  in  the  ICT  calculations. 

Also,  1  have  occasionally  made  spot  checks  of  BLAKE  runs  against  the 
new  TIGER  code;  here,  too,  the  results  are  in  excellent  agreement. 

It  is  concluded  that  BLAKE  gives  consistently  accurate  thermodynamic 
results  for  the  ideal  gas  case. 

The  non-ideal  gas  case  is  more  difficult  to  test  because  there  is 
no  generally-accepted  standard  code  that  can  be  used  as  a  standard.  The 
best  that  one  can  do  is  to  compare  one  code  against  the  other.  The  results 
are  in  Table  4.  The  agreement  is  surprisingly  good,  especially  when 
one  considers  the  difference  in  approach  between  them.  The  maximum  dis¬ 
crepancy  in  pressure  does  not  exceed  3,5%,  and  that  occurs  at  a  loading 
density  of  0.6  g/'em^,  a  figure  that  approaches  the  density  of  some 
liquids!  The  column  labeled  gives  the  ratio  of  the  computed  pressure 
to  the  ideal -gas  pressure  for  the  same  conditions.  At  a  loading  density 
of  0.6,  <5  is  2.1;  that  is,  the  correction  factor  for  the  ideal  pressure 
is  110%. 

13 

Klein  et  al.  have  proposed  a  new  equation  of  state  specifically 
for  ballistics.  Their  treatment  includes  a  term  for  the  interaction 
between  water  and  carbon  dioxide,  a  specific  interaction  that  is  ignored 
in  BLAKE.  It  will  be  interesting  to  see  the  magnitude  of  the  change  that 
this  equation  will  produce  in  ballistic  calculations;  unfortunately,  in¬ 
serting  the  equation  into  BLAKE  has  turned  out  to  be  far  from  a  routine 
matter*4. 


12 

D.  R.  Stull  and  H.  Prophet ,  Eds .,  "JANAF  Thermochemical  Tables," 

2nd  ed. ,  NSRBS-NBS  37  (June,  1971);  also  looseleaf  revisions  and  additiems 
from  the  Thermal  Research  Laboratory ,  The  Dow  Chemical  Company ,  Midland , 

MI  (M.  Chase,  Director). 

23 

E.  G.  Powell,  G .  Wilmot,  L.  Hoar,  and  M.  Klein,  "Equations  of 
State  and  Thermodynamic  Data  for  Interior  Ballistics  Calculations, " 
in  H.  Krder  and  M .  Surmerfield,  Interior  Ballistics  of  Guns,  Progress 
in  Astronautics  and  Aeronautics,  vol.  66,  AIAA ,  New  York  (1979). 

Klein,  (National  Bureau  of  Standards) ,  private  conrmmication  to 
E.  Freedman:  February  1980 . 
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TABLE  3. 


Name 

COMPARISON 

Temperature  (K) 

BETWEEN  ICT  CODE  5 
IDEAL  GAS  CASE 

Pressure  (MPa) 

BLAKE 

Impetus 

(J/g) 

ICT 

BLAKE 

ICT 

BLAKE 

ICT 

BLAKE 

One 

0.2 

2265 

2270 

173.6 

174.0 

867.8 

870.0 

0.4 

2284 

2290 

347.6 

348.5 

869.0 

871.3 

0.6 

2306 

2311 

522.2 

523.7 

870.3 

872.8 

Two 

0.2 

2558 

2560 

197.4 

197.5 

987.2 

987.6 

0.4 

2563 

2564 

394.6 

394.9 

986.6 

987.1 

0.6 

2568 

2568 

591.6 

592.0 

986.0 

986.5 

Three 

0.2 

3220 

3223 

221.8 

222.0 

1109. 

1110. 

0.4 

3229 

3252 

444.2 

444.7 

1111. 

1112. 

0.6 

3233 

3236 

666.7 

667.6 

1111. 

1113. 

Four 

0.2 

3816 

3814 

235.7 

235.7 

1178 

1178. 

0.4 

3865 

3865 

475.5 

475.7 

1189 

1189. 

0.6 

3890 

3891 

716.4 

716.9 

1194 

1195. 

Five 

0.2 

2600 

2609 

192.4 

193.1 

962.0 

965.5 

0.4 

2604 

2613 

384.6 

386.1 

961.6 

965.3 

0.6 

2608 

2617 

576.7 

579.0 

961.2 

965.1 

Pending  this  development,  one  can  take  some  comfort  from  the  agreement 
between  ELAKF.  and  the  ICT  code;  it  would  be  indeed  surprising  if  both  codes 
were  wildly  wrong.  In  my  opinion,  the  agreement  mostly  shows  the  pro¬ 
fundity  of  Hirschfelder' s  insight®  of  years  ago:  the  computations  are 
insensitive  to  the  differences  in  models  in  the  range  of  temperatures 
that  is  ballistically  interesting. 

B.  Experimental  Test 

The  over-all  conclusion  must  be  that  BLAKE  results  are  worthy  of  cre¬ 
dence,  but  the  desirability  of  experimental  confirmation  is  not  thereby 
lessened.  The  principal  difficulty  with  such  a  test  is  the  determination 
of  the  true  equilibrium  pressure,  the  maximum  pressure  that  would  be 
reached  in  a  closed  bomb  if  there  were  no  heat  loss.  This  heat  loss  causes 
a  drop  in  the  measured  pressure.  In  1931  Crow  and  Grimshaw^®  published 
an  important  paper  which  reported  such  an  attempt.  Their  pioneering  work 
necessarily  had  to  use  a  simplified  model  of  the  heat  loss  mechanism  in 
order  to  achieve  results.  A  decade  later  Kent  and  Vinti^  commented  that 
the  thermodynamic  calculation  of  the  equilibrium  temperature  was  theore¬ 
tically  more  sound;  they  used  this  temperature  in  a  new  model  of  heat 
loss  and  were  thus  able  to  compute  a  corrected  pressure.  All  later  workers 
have  followed  this  lead. 

In  1953  Vest*"  reasoned  that  heat  loss  in  a  closed  bomb  should  de¬ 
crease  with  surface-to-volume  ratio.  He  made  measurements  in  bombs  con¬ 
taining  added  metal  pieces  and  then  computed  a  corrected  pressure  by 
extrapolation  of  the  measured  pressures  back  to  a  surface-to-volume  ratio 
of  zero.  Vest's  compositions  are  listed  in  Appendix  A. 

Table  5  presents  the  final  pressures  computed  by  Crow  and  Grimshaw 
from  their  measurements  and  the  pressures  computed  with  BLAKE  for  their 
compositions.  These  compositions  are  also  listed  in  Appendix  A.  These 
calculations  were  made  without  considering  the  contribution  of  the  ace¬ 
tylene/air  igniter  that  was  used,  which  would  add  0.2  MPa  or  less  to  the 
results. 


15  A.  D.  Crow  and  V.  E.  Grimshaw,  "The  Equation  of  State  of  Propellant 
Gases,"  Phil.  Trans.  Hoy.  Soa ,  A  230 ,  30  (1931). 

2  6 

R.  II.  Kent  and  J.  P.  Vinti,  "Cooling  Corrections  for  Closed  Chamber 
Firings,"  BRL  R  281  (September,  1942).  (AD  492852) 

17 

D.  Vest,  "On  the  Transfer  of  Heat  in  High-Pressure  Combustion  Chambers , " 
unpublished  M.  S.  Thesis,  The  Johns  Hopkins  University,  Baltimore, 

Maryland  (1953). 
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The  agreement  at  loading  densities  above  0.07  g/cr  is  as  good  as 
one  could  ask.  Crow  and  Grimshaw  themselves  noted  that  their  treatment 
of  the  heat  loss  was  questionable  at  the  lower  densities.  Also,  Tent 
and  Vinti's  correction  is  much  larger  in  that  region;  their  correction 
is  much  smaller  at  the  higher  densities,  and  their  final  pressures  art: 
only  slightly  lower  than  Crow  and  Grimshaw’ s. 

Vest's  results  are  shown  in  Table  6;  they  are  much  less  extensive 
than  Crow  and  Grimshaw' s.  The  table  also  gives  the  BLAKE  results.  The 
agreement  is  good  only  at  the  loading  density  of  0.15.  The  poorer  agree¬ 
ment  at  the  two  lower  densities  may  indicate  a  failure  of  Vest’s  model. 

The  present  state  of  the  art  in  closed-bomb  experimentation  makes 
it  easier  to  conduct  experiments  at  even  higher  pressures  than  those 
used  by  Crow  and  Grimshaw.  Such  experiments  should  certainly  be  under¬ 
taken;  also,  the  Crow-Crimshaw  treatment  at  higher  loading  densities 
deserves  another  look. 


IV.  APPLICATION  OF  'BLAKE'  TO  THE  AUTOCAP  EXERCISE 

An  excellent  example  of  the  utility  of  BLAKE  :s  given  by  its  appli¬ 
cation  to  the  AUTOCAP  exercise.  One  of  the  purposes  of  this  exercise 
was  to  determine  the  effects  of  variation  in  composition  on  the  perfor¬ 
mance  of  propellants.  To  this  end,  many  lots  of  M-6  propellant  were 
made  with  deliberate  changes  from  the  accepted  specification.  Among 
numerous  other  properties,  the  impetuses  of  these  lots  were  determined; 
also,  the  impetus  of  a  reference  lot  of  M-6  was  measured.  The  ratio  of 
the  measured  impetus  of  a  lot  to  the  measured  impetus  of  a  standard  or 
reference  lot  is  called  the  relative  force  (RF). 

Table  7  gives  the  measured  and  computed  impetuses  and  the  corre-' 
sponding  relative  forces.  The  average  deviation  of  the  9  lots  is  0,6 
with  a  standard  deviation  of  0.7.  This  means  that  the  mean  deviation 
is  not  significantly  different  from  zero.  This  calculation  has  not  been 
attempted  with  another  code;  it  would  be  interesting  to  do  so. 

Another  interesting  and  useful  application  of  BLAKE  is  the  computa¬ 
tion  of  the  effect  on  the  impetus  of  small  changes  in  composition  or  of 
ingredient  properties.  This  quantity  can  be  called  the  parametric  sen¬ 
sitivity  to  distinguish  it  from  other  sensitivities  (e.g.,  ignition 
sensitivity)  to  which  it  is  unrelated. 

Mathematically  this  quantity  is  9P/9Q,  which  can  be  approximated 
by  AP/AQ.  Here,  P  is  a  thermodynamic  property  such  as  impetus,  maximum 
pressure,  temperature,  or  ratio  of  heat  capacities,  and  Q  is  the  amount 
of  an  ingredient  or  a  property  of  an  ingredient.  Thus  one  can  calculate 
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the  rate  of  change  of  maximum  pressure  with  respect  to  the  amount  of 
dinitrotoluene,  or  the  rate  of  change  of  the  ratio  of  heat  capacities 
with  respect  to  the  heat  of  formation  of  nitrocellulose,  etc. 


The  partial  derivative  notation  implies  that  all  other  constituents 
are  kept  constant.  This  '".nnot  actually  be  so,  since  it  is  only  percen¬ 
tage  composition  that  counts.  Hence  the  partial  derivative  must  be  com¬ 
puted  with  the  additional  constraint  that  the  sum  of  all  of  the  constit¬ 
uents  is  a  constant.  In  computing  AP/AQ,  the  effect  of  the  change  in 
Q  on  the  other  constituents  was  in  every  case  distributed  over  these 
other  materials  in  proportion  to  their  percentage  in  the  original  formu¬ 
lation. 

Since  impetus  or  maximum  pressure  cannot  be  readily  expressed  as  a 
function  of  composition  in  closed  form,  the  derivative  3P/3Q  must  be 
approximated  by  a  discrete  formula.  The  simplest  correct  one  is 

3P/3Q  =  (l/2h)  [P(Q  +  h)  -  P(Q  -  h)],  (27) 

where  h  is  a  small  increment  in  Q.  This  formula  approximates  the  deri¬ 
vative  at  the  given  value  of  Q.  In  order  to  decrease  the  amount  of  com¬ 
putation,  the  simpler  but  less  correct  formula: 

3P/3Q  «  (1/h)  [P(Q  +  h)  -  P CQ) 3  (28) 

was  used.  The  increment  h  was  generally  chosen  to  be  a  convenient  size, 
generally  10%  of  the  1  ise  value  of  Q. 

The  results  of  the  computations  are  given  in  Table  8.  The  sig¬ 
nificance  of  these  numbers  is  best  seen  by  looking  at  the  last  columns 
of  this  Table.  It  gives  the  approximate  percentage  change  in  impetus, 

F,  or  y  that  would  oe  brought  about  by  a  one  percent  change  in  that  ingredient. 
In  the  case  of  the  nitrogen  content,  the  number  is  the  percentage  change 
in  F  for  a  change  of  0.01  in  the  percentage  of  N. 


Increases  in  either  the  amount  of  NC  or  the  fraction  of  nitrogen 
in  the  NC  cause  an  increase  in  the  impetus;  addition  of  all  the  other 
components  causes  a  decrease.  The  reason  for  this  is  clear:  increasing 
either  NC  or  tho  fraction  of  nitrogen  increases  the  energy  of  the  pro¬ 
pellant,  hence  the  positive  sign.  Addition  of  the  other  substances 
caused  an  effective  decrease  in  the  amount  of  NC,  and  hence  decreases 
the  impetus.  Note  that  the  magnitudes  of  almost  all  of  these  quantities 
are  not  very  different  from  each  other.  This  also  is  a  result  of  the 
fact  that  the  principal  effect  of  adding  substance  Q  is  to  decrease  *he 
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BLAKE  COMPUTATIONS  ON  CROW-GRIMSHAW'S  PROPELLANTS 
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TABLE  6 


VEST’S  EXPERIMENTS 


Composition 

No.  1 

P 

Vest 

(MPa) 

BLAKE 

=  0.05 

43.7 

48.9 

0.10 

96.8 

104.2 

0.15 

159.4 

166.7 

Composition 

No.  2 

-  0.05 

54.2 

57.7 

.10 

117.9 

122.7 

.15 

193.2 

19S.0 

2S 


while  adding  it  decreases  the  percentage  of  NC,  the  net  decrease  in 
potential  energy  content  is  less  because  of  the  contribution  that  DNT 
makes. 


Knowing  the  parametric  sensitivities  is  not  sufficient.  We  also 
need  to  know  how  closely  the  level  of  a  particular  ingredient  can  be 
controlled.  A  measure  of  this  control  is  the  standard  deviation  of 
each  ingredient  in  successive  lots.  Table  9  gives  such  data*®. 


Together  the  numbers  in  Tables  8  and  9  can  give  an  estimate  of  the 
variance  of  the  muzzle  velocity.  We  consider  an  "ideal"  gun,  one  that 
has  no  energy  loss  due  to  heat  conduction,  engraving  forces,  etc.,  and 
one  in  which  the  propellant  is  completely  burnt  at  shot  start.  For  such 
a  gun,  the  muzzle  velocity,  V,  is 


V“  =  (2cF/W(y 


l))d  -  rY_1) 


where  r  is  the  ratio  of  the  volume  of  the  chamber  to  the  volume  of  the 
gun,  c/w  the  ratio  of  charge  weight  to  projectile  weight,  and  F  is  the 
impetus.  Variations  in  the  components,  Qj,  cause  variations  in  F  and 
Y,  which  in  turn  cause  variations  in  V;  that  is, 

V  =  V [ F (Qi ) ,  y(Q±)] ,  (30) 

where  stands  for  the  amounts  of  each  of  the  components.  Let  V°  be 
the  average  values  of  V.  Then  the  effect  of  deviations  from  the  averages 
of  the  Qi  can  be  expressed  by  a  Taylor  series: 


V  =  V°  ♦  T  —  +  —  ^  (Q  -  Q  •) 

1  3F  3Q.  3y  3Q.  lVi  *i 


The  variance  of  V  is  defined  bv  the  expectation  of  the  quantity 
(V  -  V°)2, 


sv2  =  e(V  -  V°)2 


This  quantity  will  now  be  calculated. 


i.  W.  May,  (Ballistic  Research  Laboratory) ,  private  communication 
to  E.  Freedman  (1975). 
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TABLE  7 


RELATIVE  FORCE  OF  AUTOCAP  PROPELLANTS 


Lot  RF  RF  Deviation 

_ Measured _ Computed _ 


A 

100.98 

100.25 

-0.73 

B 

102.03 

100.34 

-1.69 

C 

100.72 

99.72 

-1.00 

D 

102.05 

101.01 

-1.04 

E 

99.78 

99.40 

-0.38 

F 

101.23 

101.98 

0.75 

G 

95.26 

94.43 

-0.85 

H 

95.08 

94.80 

-0.28 

67994 

100.54 

99.29 

-1.25 

65306* 

100.00 

100.00 

— 

★ 

Reference  Lot 

Mean  deviation  =  -0.6 
Standard  deviation  of 

mean  =  0.7 
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TABLE  8 


Variable 

NC 

DBP 


DNT 

%  N  in  NC 


PARAMETRIC 

SENSITIVITIES  OF 

AUTOCAP  LOTS 

3F 

3Q 

3Q 

1  3F 

1  3Y 

Y -1  3Q 

3960 

-0.0012 

+0.013 

-0.0047 

-6150 

-0.020 

M) 

-6970 

-\,0 

-0.022 

'V.O 

-6380 

'uO 

-0.021 

^0 

-4400 

'ViO 

-0.014 

-vO 

-5470 

^0 

-0.018 

>v0 

-1220 

+0.0018 

-0.0039 

+0.0071 

+8100 

AjQ 

+0.026 

%o 
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TABLE  9 


VARIABILITY  OF  INGREDIENTS 
IN  AUTOCAP  M-6 


Abbreviation 

Substance 

Average 

Sigma 

%  N 

Percent  N 

13.146 

0.0114 

NC 

Nitrocellulose 

84.654 

0.4936 

DNT 

Dinitrotoluene 

10.137 

0.4930 

DBP 

Dibutyl  Pht'nalate 

5.208 

0.1478 

DPA 

Diphenyl  amine 

1.009 

0.0644 

k2so4 

Potassium  Sulfate 

1.086 

0.1161 

h20 

Water 

0.553 

0.0931 

ETOH 

Alcohol 

0.365 

2.2466 

+ 


(33) 


i>j 


2 

s 

n 
=  V 

(*3V 

3F 

3V 

+  — 

3Y 

V 

l 

i 

L9f 

3Q. 

3Y 

3Q. 

’3V  3F 

3  V 

3y 

3V 

If.  + 

3V 

3F  3Q.+ 

3y 

3Q. 

3F 

3Y 

eCQi 


Q±°) 


2  2 
The  quantity  e(Qj  -  Qj°)  is  the  variance  of  Qj,  denoted  s  (Q^). 

The  quantity  cCQ^  -  Q^°)  (Q-  -  Qj°)  is  the  covariance  of  and  Q  j , 

denoted  by  cov(Q^,  Qj). 


The  amounts  of  the  several  ingredients  are  subject  to  random  varia¬ 
tions;  hence  s^(Q^)  j£  0.  But  there  is  no  special  reason  to  believe  that 
there  is  any  relation  between  Qi  and  Qj ;  it  is  just  as  likely  that  one 
of  them  will  decrease  while  the  other  increases  as  it  is  that  both  will 
change  in  the  same  direction.  Hence  it  is  a  reasonable  assumption  to 
set  cov(Qi,  Qj)  =0. 


Differentiation  gives 


3V/3F  =  kf1/[2V(Y  -  1)],  and 
3V/3y  =  -kFf2/2V(Y  -  l)2], 

V  —  1 

where  f^  =  1  -  r'  , 

f2  =  fi  +  (Y  ~  1)  Un  r)rY_1, 
and  k  =  (2gc/W)1/2 

There  results  after  some  algebra. 


(34) 

(35) 

(36) 

(37) 

(38) 


sv2  =  (V/2fx)2 


l|W 


3F  ,  2 


3Qi 


f2V 

2  y- 


M.  1 2* 

3Qi  - 


s2(Qi) 


(39) 


A  numerical  example  is  enlightening.  Choose 
c  =  15  kg, 

W  a  45  kg, 

r  =  0.2,  y  =  1.253,  and 
F  =  106  J/kg. 
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Then  V  =  940  m/s,  and  sy  =  19  m/s,  or  about  2%.  This  is  a  significant 
variation.  Examination  of  the  numbers  that  go  into  sy  shows  that  by  far 
the  major  contributor  in  this  case  is  the  variability  in  the  alcohol. 

This  result  may  appear  paradoxical  since  the  alcohol  contributes 
virtually  nothing  to  the  impetus.  As  Sherlock  Holmes  might  have  said, 
that  is  just  the  point.  The  alcohol  contributes  nothing;  an  increase  in 
its  amount  must  come  at  the  expense  of  the  nitrocellulose,  which  does 
contribute  a  great  deal.  If  the  variability  of  the  alcohol  had  been  only 
0.5  (the  same  as  NC),  the  resultant  sy  would  have  been  4  m/s,  or  about 
0.4%,  a  tolerable  amount. 


V.  SOURCES  OF  DATA 

BLAKE  is  presently  dimensioned  to  permit  the  simultaneous  considera¬ 
tion  of  29  gases  and  10  condensed  phases  (liquids  or  solids  or  both)  at 
one  time;  these  species  may  contain  up  to  40  elements.  In  practice  cer¬ 
tain  constraints  (discussed  in  Sec.  VII)  must  be  considered. 

A.  Elements 

BLAKE  presently  considers  16  elements:  C,  H,  N,  0,  B,  Al,  F,  Na, 

Ti,  He,  Ar,  Ba,  K,  S,  Cl,  and  Pb.  Users  may  add  another  24  without  having 
to  redimension  any  variables. 

B.  Product  Species 

Thermochemical  data  for  120  product  compounds  (98  gases  and  22  con¬ 
densed  species)  are  included  in  the  library  furnished  with  BLAKE.  The 
data  for  these  compounds  are  taken  without  exception  from  the  JANAF 
Tables-*^.  They  are  stored  in  the  form  of  the  coefficients  of  a  linear 
combination  of  functions  of  T/1000  (T  =  absolute  temperature);  these  co¬ 
efficients  were  determined  by  a  least-squares  fitting  and  reproduce  the 
data  over  the  range  400  K  -  5000  K.  (See  also  Sec.  VII  C  1.) 

The  Lennard-Jones  parameters  were  taken  from  the  compilation  of  Reid 
and  Sherwood*®,  In  the  many  cases  when  no  data  are  available,  the  program 
automatically  substitutes  e  =  300  K  and  a  =  3.5  x  10“*®  meter.  The  user 
can  change  these  default  values  permanently  at  compile  time  by  changing 
them  in  Subroutine  LJPAR;  or  temporarily  at  run  time  by  using  the  instruc¬ 
tion  LJP  (Sec.  VII  C  2). 


heidj  H.  C.  and  Sherwood ,  T. K. ,  "The  Properties  of  Gases  and  Liquids 
2nd  ed.y  McGraw-Hill ,  New  York  (1966). 
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C.  Ingredients 


BLAKE  contains  enthalpies  of  formation  for  349  different  species  in 
a  DATA  statement.  Of  these,  300  are  nitrocelluloses  with  percentages  of 
nitration  varying  from  11.00  to  14.00.  These  data  are  based  on  a  re¬ 
working  of  those  published  by  Jessup  and  Prosen2^  (including  the  correc¬ 
tion  of  an  obvious  misprint).  A  least  squares  fitting  of  these  data  gives 

Ah  =  -1406.22  +  6231.86  f  ,  (40) 


where 


Ah  =  enthalpy  of  formation  (cal/g). 


and 


f  =  weight-fraction  of  nitrogen  in  the  nitrocellulose. 


Nitrocellulose  with  100  f  weight -percent  nitrogen  has  the  empirical 
formula  X05+2X^X»  where  x  is  the  number  of  nitrogen  atoms  in  a 

monomeric  unit,  given  by 


x  =  162.1430  f/ (14. 0067  -  44.9975  f) 

The  molecular  weight  corresponding  to  this  empirical  formula  is 


(41) 


MW  =  44.9975  +  162J430f  (42) 

The  combination  of  these  formulas  with  Jessup  and  Prosen's  data  gives  the 
enthalpies  of  formation  of  all  nitrocelluloses  as  a  function  of  f. 

The  other  50  compounds  are  49  of  the  most  common  ingredients  used  in 
the  formulation  of  solid  and  liquid  propellants.  (One  compound  appears 
twice.)  The  names  of  these  compounds  and  the  abbreviations  or  acronyms 
adopted  for  them  are  listed  in  Table  10. 

As  much  as  possible,  the  acronyms  or  abbreviations  were  chosen  to 
agree  with  customary  and  traditional  usage.  Ethanol  appears  twice,  one 
as  ALC  for  those  who  prefer  the  common  abbreviation,  and  once  as  ETOH  for 
those  who  appreciate  the  subtleties  of  organic  nomenclature. 


20 

R,  S.  Jessup  and  E.  J.  Prosen ,  " Heats  of  Combustion  and  Formation 
of  Cellulose  and  Nitrocellulose  ( Cellulose  Nitrate) J.  Res.  National 
S  " eau  of  Stds.  44 ,  387-393  (1950). 
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TABLE  10 


Chemical  Name 


Acetone  (dimethyl  ketone) 

Alcohol  (ethanol) 

Amyl  phthalate 
Akardite  2 
Ammonium  nitrate 

Barium  oxide 
Barium  nitrate 
Basic  lead  oxide 

Carbon  (graphite) 

Cellulose  diacetate 

Cellulose  nitrate  —  see  nitrocellulose 
Cellulose  triacetate 
Cryolite  (sodium  aluminum  fluoride) 
Cyclotetramethylene  tetranitramine 
(octogen) 

Cyc lo- 1,3, 5 -trimethyl ene-2 , 4 , 6- trini tramine 
(hexogen) 

Decalin 

Dibutyl  phthalate 
Diethyleneglycol  dinitrate 
2,2' -dinitrodiphenylamine 
Dinitrotoluene 
Diphenyl amine 

Ethanol 

Ethanolamine  nitrate 
Ether 

Ethyl  centralite 
Ethylenediamine  dinitrate 
Ethylene  glycol 
Ethyl  ether 

Glycerin 

Graphite 

H2O  (water) 

Hexogen 

HMX 

Hydrazine 
Hydrazine  nitrate 
Hydroxylaramonium  nitrate 


Abbreviation 

ACETON 

ALC 

AMYLPH 

AKAR2 

AN 

BAO 

BANITR 

PB2C04 

C 

CEDAC 

CETAC 

CRY 

HMX 

RDX 


DECA 

DBP 

DEGDN 

NNDP 

DNT 

DPA 

ETOH 

EOAN 

ETHER 

EC 

EDDN 

EGLY 

ETHER 

GLY 

C 

H20 

RDX 

HMX 

N2H4 

HN 

HAN 
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Chemical  Name 


Abbreviation 


Isopropylammonium  nitrate  IPAN 

Lead  oxide  (see  also  basic  lead  c:t;.de)  PBO 

Metriol  METRIC 

Methyl  central ite  MC 

Nitrocellulose  with  vw.xy  NCvwxy 

percent  nitrogen 

2-nitrodiphenylamine  NDPA 

Nitric  acid  HN03 

Nitrogen  N2 

Nitroglycerin  N6 

Nitroguanidine  NQ 

Octogen  HMX 

Oxygen  02 

Pentaerythritol  pentanitrate  PETN 

Potassium  nitrate  KN 

Potassium  sulfate  KS 

Potassium  tetrafluoroborate  Ki.’F4 

RDX  RnX 

Sodium  aluminum  fluoride  CRY 

Triacetin  TRIAC 

Triaminotrinitrobenzene  TATB 

Triethyleneglycol  dinitrate  TECjJN 

Trimethylammonium  nitrate  TMAN 

Unsymmetrical  dimethyl  hydrazine  UDMH 

Water  H20 
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The  formulas  and  enthalpies  of  formation  of  these  ingredients  are 
listed  in  Appendix  B.  Wherever  possible  they  were  taken  from  the  author¬ 
itative  compilation  of  Cox  and  Pilcher^*.  Data  for  the  inorganic  com¬ 
pounds  were  taken  either  from  the  NBS  compilation^?  or  from  Lange's 

Handbook^ . 


VI.  IMPLEMENTING  'BLAKE' 

A.  Machine-Dependent  Features 

BLAKE  is  written  almost  exclusively  in  the  extended  set  of  American 
National  Standard  FORTRAN.  There  are  six  functions  and  one  feature  in 
it,  however,  that  are  system-  or  computer-dependent.  These  functions 
(CLOCKS,  SECOND,  DAY,  ENCODE,  DECODE,  .AND  EOF)  are  used  only  in  convenience 
features  and  can  be  readily  omitted.  The  machine-dependent  feature,  word 
length,  may  cause  a  little  more  work. 

1 .  Clock. 

It  is  desirable  but  not  essential  that  the  program  be  capable  of 
obtaining  the  current  time  from  the  system  clock.  As  furnished,  the  pro¬ 
gram  contains  Subroutine  CLOCKS(T).  This  subroutine  returns  the  time  in 
seconds  from  the  previous  midnight.  BLAKE  obtains  this  time  from  a  system 
function  named  SECOND (T). 

2 .  Date. 

It  is  desirable  but  not  essential  that  the  program  be  capable  of 
obtaining  the  current  date  from  the  system  clock  or  calendar.  As  furnished, 
the  program  contains  Subroutine  DAY(IDAY).  IDAY  is  an  array  of  20  alpha¬ 
numeric  characters  that  is  later  printed  using  a  format  of  20A1.  BLAKE 
inserts  the  current  date  in  the  first  12  places,  and  blanks  in  the  other 
8.  It  does  this  using  the  machine-dependent  functions  ENCODE  and  DECODE. 

If  you  want  this  feature,  you  will  have  to  make  appropriate  implementations 
of  these  functions. 


21 

J.  D.  Cox  and  G.  Pilcher,  " Thermochemistry  of  Organic  and  Organo- 
metallic  Compounds ,"  Academic  Press,  London  (1970). 

22 

D.  D.  Wagman,  et  al.,  " Selected  Values  of  Chemical  Thermodynamic  Proper¬ 
ties,"  NBS  Tecbr.ical  Note  270-3,  Washington,  DC  (1968). 
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J.  A.  Dean,  ed. ,  "Lange’s  Handbook  of  Chemistry ,"  12th  ed. ,  McGraw- 
Hill,  New  York  (1979). 
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3.  EOF. 


It  is  convenient  but  not  necessary  to  check  for  the  occurrence 
of  an  end-of-file  mark  in  the  input.  BLAKE  does  this  in  Subroutine  CARDRD 
after  each  READ  instruction  by  use  of  the  CDC  function,  EOF(LI),  where 
LI  is  the  logical  unit  number  of  the  input  device.  You  will  probably  have 
to  change  it  for  your  particular  computer. 

The  standard  FORTRAN-77  statement 

READ  (unit,  format,  END=435) 

is  already  in  the  program  as  a  comment  card.  If  your  computer  will  rec¬ 
ognize  it,  by  all  means  use  it. 

4.  Word  Length. 

TIGER,  and  hence  BLAKE,  proceed  by  interpreting  a  series  of  in¬ 
structions.  The  detailed  form  of  these  instructions  is  given  in  Section 
VI;  here  it  is  only  necessary  to  note  that  each  separate  field  in  an  in¬ 
struction  is  converted  by  the  program  (in  Subroutine  CARDRD)  to  a  real 
number  containing  from  2  to  12  digits  (2  digits  per  character) .  For 
example,  the  symbols 

P,  H20,  17E-3,  ana  K2S04$ 

become 

30.,  220229.,  0107191303.,  and  250233290447., 

respectively.  Some  of  thc^e  numbers  are  used  as  such;  others,  such  as 
220229  ('H20’)  are  eventually  converted  back  to  their  original  alpha¬ 
numeric  form  for  printing.  The  conversion  routine.  Subroutine  CONVRT, 
starts  by  converting  the  real  number  into  integer  form.  Depending  on 
the  number  of  characters  permitted,  the  routine  may  have  to  handle  from 
6  to  12  digits  without  the  loss  of  even  the  least  significant  digit. 

CDC  computers  can  handle  12  digits  without  any  loss,  so  BLAKE 
is  set  to  accept  6  characters  in  Subroutine  CARDRD.  Most  other  computers 
cannot  handle  more  than  8  digits  (4  characters)  and  hence  require  changes 
in  CARDRD. 

The  change  is  simple.  Line  100  currently  reads 
IF  (IFLAG  -  6)  270,  280,  290 

All  that  is  needed  is  to  change  the  '6'  to  a  '4'.  If  your  computer  can 
only  handle  6  or  7  digits,  the  '6'  becomes  a  * 3 * . 
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Also,  two  changes  are  required  in  Subroutine  LIBRAY 
immediately  after  statement  1010  which  now  reads 


First 


the  line 


IF  (STR(J,2)  - 
must  be  changed  to 

IF  (STR(J, 2)  - 
Secondly,  the  line  immediately 
IF  (STC(J,2)  - 
must  be  changed  to 

IF  (STC(J,2)  - 

Two  changes  are  needed  in 
the  number 


3400000000.)  1020,  1020,  1060 

30000000.)  1060,  1060,  1020 
after  statement  1120 
3400000000.)  1140,  1140,  1130 

30000000.)  1130,  1130,  1140 

Subroutine  BLAKEM.  Statement  390  contains 


it  must  be  changed  to 


2318201526.  ; 


231820. 


Statement  420  contains  the  number 

362332231526.  ; 


it  must  be  changed  to 


362332. 

Finally,  changes  will  be  needed  in  the  library  itself. 

1.  The  words  SOLID,  LIQUID,  AND  CONDEN  must  to  changed  to  SOLI, 
LIQU,  and  COND,  respectively. 

2.  Many  of  the  constituents  in  the  library  have  six-character 
names;  these  must  all  be  changed  to  four  characters  each.  Be  careful  when 
making  the  changes  to  note  two  things: 

a.  All  names  must  be  unique;  it  is  regrettably  easy  to  in¬ 
troduce  duplicates  when  shortening  the  existing  names. 

b.  If  possible,  make  the  changes  with  a  computer  editor. 
Some  of  the  names  may  occur  14  times;  it  is  so  very  easy  to  miss  one 
manually . 
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B. 


User-Dependent  Aspects 


There  are  some  things  in  BLAKE  that  are  matters  of  user-  or  installa¬ 
tion-preference,  such  as  the  designation  of  input  and  output  units, 

1.  The  input  and  output  units  for  BLAKE  are  named  symbolically 

LI  and  L0  everywhere,  with  three  exceptions.  LI  and  L0  are  set  equal  to 
5  and  6,  respectively,  in  the  main  program.  The  three  exceptions  occur 
in  Subroutines  DVDINT,  MATINV,  and  TICTIC,  where  the  WRITE  unit  is  ex¬ 
plicitly  named  unit  6.  You  may  need  to  change  these  designations. 

2.  The  program  forms  and  then  uses  a  binary  file  on  unit  LIB,  which 
is  set  to  7  in  the  main  program.  Again,  you  may  prefer  a  different 
number . 

3.  The  units  for  input  quantities  in  BLAKE  are: 

energy:  calories; 
temperature:  Kelvins; 
density:  grams  per  cubic  centimeter; 
pressure:  atmospheres. 

These  units  can  be  readily  changed  by  altering  one  or  more  of  the  con¬ 
stants  PCON,  VCONj  TMCON,  TACON,  and  HCON  in  the  main  program.  These 
are  used  throughout  the  program  as  conversion  factors  for  pressure, 
volume,  temperature  (two  factors),  and  energy,  respectively.  They  are 
now  set  to  1.0,  1.0,  0.0,  1.0,  and  1,0. 

For  example,  if  you  wanted  to  change  input  units  from  calories 
and  Kelvins  to  joules  and  degrees  Fahrenheit,  the  following  changes 
would  be  made  before  compiling  the  program: 

HCON  =  4.184, 

TbCON  =  0.555556,  and  TAC0H  =  255.375. 


VII.  USING  'SLAKE' 

A.  Job  Control  Cards  and  the  First  Two  Runs 


All  job  control  cards,  system  cards,  or  t.heir  equivalent,  are  the 
responsibility  of  the  user. 

In  what  follows,  it  will  be  assumed  that  the  program  has  been  com¬ 
piled  and  goaded  without  error.  The  very  first  time  the  program  is 
executed,  it  is  necessary  to  form  the  binary  library  that  will  there¬ 
after  be  used  in  every  run.  Appropriate  job  control  instructions 
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must  be  supplied  to  assign  a  file  for  this  library,  which  the  program 
assumes  will  be  on  unit  7,  although  you  can  readily  change  this  assign¬ 
ment.  (See  Section  VI. B) 

The  BLAKE  library,  which  is  supplied  along  with  the  program,  is  the 
input  for  this  first  run.  The  output  produced  consists  of  the  binary 
library  on  unit  7,  and  a  human-readable  output  on  the  chosen  output  device 
(unit  6  unless  you  have  previously  changed  it) .  This  latter  output  should 
be  carefully  scrutinized  for  error  messages,  which  will  be  the  first  part 
of  the  output.  Assuming  that  all  of  the  necessary  changes  described  in 
Section  VI  were  made  prior  to  this  run,  there  should  not  be  any  error 
messages  from  the  operating  system  at  this  point.  There  may,  however, 
be  error  messages  from  the  program  calling  attention  to  problems  with 
the  library.  The  most  common  one  is  TWO  OR  MORE  CONSTITUENTS  HAVE  THE 
SAME  NAME,  which  will  occur  if  the  reduction  of  the  names  of  products 
to  four  or  fewer  characters  was  not  done  carefully. 

Following  the  error  messages  (if  any)  will  be  a  complete  listing  of 
the  library  in  human-readable  form.  It  should  be  painstakingly  read 
to  guard  against  unexpected  problems  (garbaged  names,  for  example) . 

If  it  looks  good,  then  you  are  ready  for  the  second  run,  which  will 
consist  of  the  test  cases  also  supplied  with  the  program.  Before  run¬ 
ning  them,  it  is  necessary  to  insert  control  cards  that  will  assign  the 
binary  library  formed  in  the  previous  run  to  unit  7. 

The  output  from  the  test  cases  should  be  compared  with  the  listing 
given  in  Appendix  C  of  this  report.  Further  discussion  of  the  output 
is  given  in  Section  VIII. 

B.  Comments  on  the  Computation  of  Equilibrium  by  TIGER 

In  order  to  shorten  the  discussion  required  for  some  of  the  instruc¬ 
tions,  a  few  comments  are  required  on  the  philosophy  of  TIGER'S  calcula¬ 
tion  of  the  equilibrium  state.  The  discussion  will  be  entirely  qualita¬ 
tive;  the  complete  technical  details  are  in  the  TIGER  manuals-^. 

TIGER  proceeds  by  considering  at  first  only  the  gaseous  constituents 
that  are  chemically-possibie  nroducts  of  the  specified  mixture.  Conden¬ 
sed  phases  are  usually  (not  uivcys)  considered  later  if  necessary.  At 
each  step  the  program  tries  to  work  with  a  set  of  N  linearly-independent 
constituents  that  are  present  in  the  largest  amounts,  where  N  is  the 
number  of  different  elements  present  in  the  mixture.  These  N  constituents 
are  called  the  'components’.  Most  of  the  computing  time  is  spent  looking 
for  the  correct  set  of  components,  which  may  change  with  conditions, (e.g. , 
pressure  or  temperature)  even  while  computing  a  single  composition. 
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In  theory  but  not  in  practice  the  program  could  permute  all  of  the 
constituents  present  until  it  finds  the  components.  The  first  attempt 
at  finding  the  components  consists  of  selecting  the  first  linearly-inde- 
pendent  set  of  constituents  as  they  are  retrieved  from  the  library  tape. 
This  order  was  determined  at  the  time  the  library  tape  was  created  and 
was  based  on  the  order  of  the  CONstituent  cards.  For  the  elements  C, 

H,  0,  and  N,  the  components  usually  turn  out  to  be  CCb,  FbO,  NS,  and  CO, 
which  are  therefore  the  first  four  constituents  in  the  BLAKE  library. 

See  also  the  discussion  of  the  instructions  QRDer  and  ORC  in  Section  VII 
C  2. 


If  the  first  selection  of  components  proves  to  be  wrong,  the  program 
does  some  permuting  among  the  constitutents.  Since  one  knows  chemically 
that  some  of  the  constituents  will  rarely  if  ever  be  present  ir.  large 
enough  concentrations  to  be  components,  the  program  makes  provision  for 
skipping  them.  Constituents  that  are  never  to  be  considered  as  components 
have  an  asterisk  (*)  punched  as  the  last  field  in  their  CONstituent  in¬ 
struction  (See  Section  VII  Cl).  When  the  program  prints  the  list  of  con¬ 
stituents  selected,  those  that  cannot  be  components  are  all  printed  after 
those  that  will  be.  The  dividing  line  between  these  two  groups  is  called 
the  'floor'.  You  have  the  capability  of  raising  or  lowering  the  floor; 
see  FLOor  in  Section  VII  C  2). 

BLAKE  is  presently  limited  to  considering  a  maximum  of  29  gaseous 
constituents  for  any  one  composition.  Error  message  No.  31  is  printed 
if  more  than  29  gaseous  constitutents  are  selected,  and  the  case  is  re¬ 
jected. 

The  library  contains  thermodynamic  d3ta  for  a  wide  range  of  species. 
Consequently,  compositions  with  five  or  more  elements  can  readily  have 
mere  than  29  constituents  chosen  by  the  program  and  thus  will  not  run. 

This  potentially  embarrassing  problem  is  solved  by  the  REJect  instruction. 
(See  Section  VII  C  2). 

Convergence  is  materially  aided  if  you  have  advance  knowledge  of  the 
major  constituents  in  the  equilibrium  composition.  Such  knowledge  is 
often  obtained  by  running  the  same  composition  with  another  code.  The 
NASA-Lewis  equilibrium  code  is  recommended  for  this  purpose.  The  NOTS 
code  used  at  the  Naval  Weapons  Center  also  works  well. 

The  NASA-Lewis  code  has  the  unsurpassed  advantage  that  for  a  given 
composition,  it  can  consider  all  of  the  chemically-pcssibic  species  con¬ 
tained  in  the  JANAF  Tables.  Thus,  one  can  run  an  obstreperous  formula¬ 
tion  (cne  containing  more  than  six  elements)  with  this  code  first,  deter¬ 
mine  the  significant  constituents,  and  then  use  th^s  information  to 
determine  REJect,  ORDer,  and  ORC  instructions  for  a  subsequent  BLAKE  run. 
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C.  General  Remarks  on  the  Instrir  Mon  Set 


BLAKE,  like  TIGER,  proceeds  by  executing  one  by  one  a  set  of  instruc¬ 
tions;  the  BLAKE  set  contains  37,  of  which  a  third  are  truly  essential, 
a  third  serve  to  change  default  values  of  diverse  parameters,  and  a  third 
are  frankly  frills,  intended  to  modify,  expand,  contract,  label,  title, 
or  comment  on,  the  output. 

To  a  considerable  extent  the  actual  order  in  which  the  instructions 
are  entered  by  the  user  is  immaterial.  Logic,  to  be  sure,  demands  that 
some  orders  precede  others:  a  heat  of  explosion  cannot  be  computed  until 
a  composition  has  been  specified;  but  a  title,  for  example,  can  be  entered 
late  or  not  at  all,  as  you  please. 

The  language  of  the  descriptions  of  these  instructions  generally 
assumes  that  they  are  created  by  being  typed  into  an  input  file  from 
a  terminal,  but  they  can  just  as  well  be  entered  from  punched  cards. 

Each  instruction  consists  of  a  key  word  or  abbreviation  followed 
by  a  variable  number  of  fields  containing  keywords  or  numbers.  Each 
field  is  set  off  from  its  neighbors  by  commas.  The  key  word  or  abbrev¬ 
iation  for  the  instruction  itself  must  appear  in  columns  1-3  of  the 
instruction;  otherwise  the  fields  are  arbitrary  in  both  length  and  place¬ 
ment  in  the  line. 

When  numbers  are  required  in  an  instruction,  they  may  be  typed  in 
free  form.  Thus,  69,  69.,  69.0,  0.69E2,  and  0.0069E04  all  represent 
the  same  number  in  an  instruction.  The  program  decides  eventually  whether 
it  needs  a  real  number  or  an  integer  and  takes  appropriate  action. 

The  general  form  of  an  instruction  is: 

INStruction,  PARI,  nl,  PAR2,  n2,  PAR3,  n3,  .  .  . ,  (PARm,  nm) 

Names  or  letters  in  CAPITALS  must  be  typed  just  that  way;  no  variation  is 
permitted.  Lower-case  letters  are  added  to  some  of  the  instruction  names 
for  mnemonic  reasons,  but  are  not  part  of  the  instructions  name  and  may 
be  omitted.  Thus,  FOR,  FORmula,  and  FORmulas  are  all  interpreted  as  the 
same  instruction  and  all  are  valid.  The  comma  marks  the  end  of  this 
field  and  separates  succeeding  fields  from  each  other.  In  a  few  cases, 
double  commas  with  at  most  blanks  between  them  have  a  special  significance; 
their  use  is  explained  in  the  applicable  instruction. 

The  quantities  PARI,  PAR2,  PAR3,  ...  ,  represent  key  words  or  letters 
and  must  be  typed  just  as  given  for  each  particular  instruction.  The 
quantities  'nl',  'n2',  ’nS1,  ...  ,  represent  numbers  that  are  associated 
with  the  keywords  or  parameters  and  may  be  typed  as  desired  or  needed. 
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Sometimes  the  parameters  must  appear  in  a  specified  order;  othertimes 
order  is  immaterial.  When  the  order  may  be  varied,  this  is  so  stated; 
if  not,  assume  that  the  items  must  appear  in  the  given  order. 


The  ellipsis,  ...  ,  indicates  that  a  variable  number  of  items  may 
be  entered.  Names  in  parentheses  are  optional,  as  noted  for  individual 
instructions. 

D.  The  Instruction  Set 

1.  Library  Instructions.  The  library  instructions  are  CONstituent, 
ELEment,  PRInt,  STC,  STG,  STR,  and  END.  These  instructions  are  used  to 
form  the  binary  library  and  are  understood  by  the  program  only  while  that 
operation  is  proceeding.  This  operation  in  turn  starts  with  the  instruc¬ 
tion,  STArt  of  library,  which  is  thus  not  a  library  instruction,  but  which 
has  no  other  purpose. 

CONstituent,  name,  GAS,  el,  nl,  e2,  n2,  e3,  n3,  ... 


CONstituent,  name,  CONDEN,  el,  nl,  e2,  n2,  e3,  n3,  ... 

This  library  instruction  defines  'name'  as  the  official  name  of 
a  constituent.  One  of  the  parameters  GAS  or  CONDEN  (for  condensed)  must 
appear.  The  quantities  ’el',  'e2',  etc.,  are  the  symbols  of  the  elements 
that  appear  in  'name',  ’nl’,  'n2',  etc.,  are  the  numbers  of  atoms  of  those 
elements  respectively.  Note  that  a  condensed  constituent  may  appear  in 
either  the  liquid  state,  the  solid  state,  or  both. 

ELEment,  sym,  atwt 

This  library  instruction  defines  •sym'  as  the  official  name  of  an 
element  whose  atomic  weight  is  'atwt'.  No  constituent  name  will  be 
accepted  by  the  program  unless  all  of  its  constituent  atoms  appear  in 
ELEment  instructions. 

Note  that  if  the  element  itself  is  to  be  allowed  as  a  constituent, 
it  must  also  appear  in  a  CONstituent  instruction. 

Examples: 

ELEment,  NA,  22.98997 
CONstituent,  NA,  GAS,  NA,  1 

If  the  latter  instruction  were  omitted,  then  sodium  could  appear  as 
NaOH  or  NaO,  etc.,  but  not  as  the  monatomic  gas,  Na. 
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END  of  Library 

This  instruction  marks  the  end  of  the  library  deck.  It  must  appear. 

As  soon  as  it  is  read,  the  program  then  proceeds  to  sort  the  library  cards, 
form  the  binary  library  tape  (or  its  equivalent  on  the  disk),  and  writes 
the  library  listing  (unless  the  latter  is  suppressed) .  Note  that  library 
instructions  may  appear  in  any  order  between  STArt  and  END;  the  program 
sorts  them  after  the  END  card  has  been  read. 

STC,  name,  SOLID,  1,  al,  a2,  a3 

STC,  name,  SOLID,  2,  a4,  a5,  a6 

STC,  name,  SOLID,  3,  a7,  a8,  a9  (and  similar  forms  for  LIQUID) 


The  STC  instruction  supplies  data  for  the  equation  of  state  of  a 
condensed  substance,  either  LIQUID  or  SOLID.  The  order  of  the  cards  is 
immaterial,  but  all  three  must  appear;  if  not,  all  of  the  numbers  are 
skipped.  (This  will  eventually  lead  to  a  run-time  error,  when  the  program 
tries  to  take  the  reciprocal  of  'al*.)  The  three  coefficients  must  appear 
in  the  correct  order  on  a  given  card,  however. 


The  coefficients  themselves  belong  to  the  empirical  equation 
V  =  Aj  +  A2P  +  a3P2, 

where  V  is  the  molar  volume  of  'name',  and  the  A^  are  given  by 

A1  =  ai  +  a2T  +  a3j2  > 

2 

A~  =  a.  +  arT  +  a,T  ,  and 
i  4  o  o 

A3  =  a?  +  agT  +  agT  . 


(45) 

(44) 

(45) 

(46) 


Here,  P  and  T  are  pressure  (atmospheres)  and  temperature  (Kelvins); 
the  a^  are  represented  by  the  'ai'  in  the  instruction.  It  is  unlikely 
that  all  9  coefficients  will  be  known  for  many  materials,  but  zeroes 

must  appear  on  the  cards  for  all  missing  coefficients. 

( 

STR,  name,  GAS,  1,  bl,  b2,  b3 

STR,  name,  GAS,  2,  b4,  b5,  b6 

STR,  name,  GAS,  3,  b7,  b8,  b9 

These  instructions  introduce  the  thermodynamic  data  for  ’name'  into 
the  library.  LIQUID  or  SOLID  may  appear  instead  of  GAS.  If  a  substance 
forms  only  a  liquid  (solid),  then  the  instructions  for  solid  (liquid) 
are  omitted  for  that  substance. 


44 


The  coefficients,  'bi',  are  related  to  the  heat  capacity  at  constant 
pressure,  C,  enthalpy,  H,  and  entropy,  S,  by  the  empirical  equations 

C  =  bj  +  b20  +  b302  +  b403  +  b5/0  +  b6/02  +  b7/©3;  (47) 

S  =  R(b1ln(0)  +  b20  +  b302/2.O  +  b4©3/3.0  -  b5/0  -  b6/2.O02  - 

67/3.O02);  and  (48) 

H  =  [RTO^  +  b_0/2.O  +  b302/3.O  +  b4©3/4.0  +  b5ln(0)/0  - 

b6/02  -  by/2.O03)  +  bg  -  AH° (298 . 15) ]/ 1000.0  (49) 

Here, 


0  =  T/1000, 

R  is  the  gas  constant,  1.98717  cal/mole-K,  and  AH°(298.15)  is  the  standard 
enthalpy  of  formation. 

These  coefficients  are  obtained  by  applying  least-squares  techniques 
to  suitably-selected  data,  preferably  those  in  the  JANAF  Tables1^. 

BLAKE  is  furnished  with  a  library  containing  data  for  98  gases,  and  22 
condensed  phases  of  16  elements.  Users  desiring  additional  data  must  do 
their  own  fittings. 

PRInt  the  Library 

This  instruction  creates  a  listing  of  the  library  on  output  unit 
LO.  The  listing  is  not  printed  unless  this  instruction  is  entered.  It 
may  appear  anywhere  between  ST Art  and  END.  The  library  listing  is  ordered 
according  to  the  order  of  appearance  of  the  CONstituent  instructions  and 
is  formatted  for  human  reference.  A  copy  should  be  made  every  time  the 
library  is  changed. 

STG,  name,  cov,  T,  t,  S,  s 

This  instruction  differs  from  the  STG  instruction  of  TIGER.  Here, 
•name*  is  the  legal  name  of  a  gaseous  constituent  as  defined  by  a  CON¬ 
stituent  instruction.  The  quantity  'cov'  is  the  BKW  covolume  of  this 
substance  (not  to  be  confused  with  the  covolume  used  in  the  Noble-Abel 
equation) . 

The  quantities  *t’  and  's'  are  the  two  Lennard-Jones  parameters, 
the  well  depth  and  the  molecular  diameter.  The  key  letters  T  and  S  may 
be  replaced  by  either  E  or  D,  respectively,  and  their  order  may  be  inter¬ 
changed.  Of  course,  ' t'  and  's'  must  be  associated  with  the  proper  key 
letter. 
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It  is  sufficient  to  supply  only  the  BKW  covolume,  but  if  one  of  the 
Lennard-Jones  parameters  is  specified,  then  so  must  the  other;  if  not, 
neither  one  is  accepted. 

2.  Computational  and  Control  Instructions.  The  instructions  are 
listed  in  alphabetic  order  by  their  three- letter  key. 

CHOose,  narael,  name2,  name3,  ... 


Example: 

CHO,  co2,  h2o,  n2,  CO,  nh3 

This  instruction  forces  the  program  to  consider  only  the  named  con¬ 
stituents  as  possible  products  of  an  equilibrium  calculation.  The  names 
listed  must  be  constituent  names  in  the  library.  The  use  of  this  instruc¬ 
tion  is  recommended  only  in  special  cases. 

CM2 

See  the  discussion  under  COMposition. 

CMT  (message) 

This  instruction  prints  the  word  COMMENT  followed  by  the  'message' 
you  have  typed,  all  on  one  line.  It  is  convenient  for  adding  notes  to 
the  output. 


Example: 

CMT  this  is  a  comment  card. 

COMposition,  namel,  al,  name2,  a2,  name3,  a3,  _  ,  (MOLE) 

This  instruction  defines  the  composition  of  a  mixture.  It  is  an 
essential  instruction,  as  all  of  the  computational  instructions  are  auto¬ 
matically  rejected  if  they  occur  before  a  COM  instruction. 

The  names  'namel',  'name2',  etc.,  must  have  previously  been  de¬ 
fined  by  FORmula  instructions;  otherwise,  the  entire  instruction  is  re¬ 
jected.  The  quantities  al,  a2,  ...,  are  the  relative  amounts  of  each 
substance.  The  order  of  the  names  is  immaterial,  but  the  correct  amount 
must  follow  each  name.  The  amounts  need  not  be  percentages;  any  basis 
is  acceptable.  Regardless  of  how  you  enter  the  composition,  however, 
the  program  will  normalize  it  and  print  it  in  both  weight-percent  and 
mole  per  cent. 

If  you  prefer  to  enter  your  composition  data  in  moles,  the  keyword 
MOLE  must  also  be  typed. 
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All  of  the  data  must  fit  into  one  80-column  line.  If  this  is  not 
feasible,  use  the  instruction  CM2  instead  of  COM.  You  may  then  extend 
the  input  data  over  two  80-column  lines.  Do  not  repeat  CM2  at  the  start 
of  the  second  line.  The  MOLE  designator  does  not  apply  to  CM2. 

Examples: 

C0M,  NC1325,  75,  NG,  25 
C0MP0S,  N2 ,  4,  02,  1,  M0LE 

DATe  (message) 

The  separator  (,)  is  not  needed  by  this  instruction.  The  date  of 
the  run  is  normally  obtained  from  the  operating  system  clock  and  printed 
on  each  output  page.  If  you  have  trouble  in  obtaining  this  information, 
or  if  you  want  to  print  your  own  date,  this  instruction  takes  the  last 
20  characters  from  'message'  and  prints  them  in  place  of  the  system  date. 
If  the  'message'  is  blank,  the  system-produced  date  is  omitted  and  nothing 
replaces  it.  Once  DATe  has  been  used,  you  cannot  return  to  using  the 
system  clock  in  the  same  run. 

The  regular  use  of  this  instruction  is  not  recommended  if  you  have 
access  to  the  system  date.  The  primary  purpose  of  the  instruction  is  to 
mark  special  days,  like  one's  birthday,  Halloween,  Friday  the  13th,  Bands¬ 
men's  Day,  etc. 

Example: 

DAT  Fourth  of  July 


DEBug,  n 

If  'n'  is  not  0,  the  debug  feature  is  turned  on;  if  zero,  it  is 
turned  off  if  it  was  on. 

The  'debug'  feature  was  taken  over  from  TIGER,  where  it  was  introduced 
for  assistance  in  developing  Subroutine  ECOMPO.  When  turned  on,  it  prints 
the  current  temperature,  pressure,  volume,  entropy,  and  energy,  the  cur¬ 
rent  mole  numbers,  and  the  current  corrections  to  the  mole  numbers  on 
each  pass  through  ECOMPO.  In  my  experience  it  is  net  presently  useful, 
and  will  be  removed  from  later  versions  of  BLAKE. 

DEStroy 

This  instruction  negates  the  effect  of  a  RETain  (which  see)  instruction. 
ECHo,  n 
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This  instruction  makes  the  program  print  out  {'echo')  each  instruction 
as  it  is  encountered.  The  feature  is  turned  on  when  'n*  is  non-zero,  and 
turned  off  when  'n'  is  0. 

EXPlosion,  V,  vol,  (EOF,  e) 

This  instruction  executes  a  constant-volume  (isochoric)  burning  cal¬ 
culation  on  a  previously-defined  mixture.  The  specific  volume,  vol,  is  in 
cubic  centimetres  per  gram  unless  you  have  changed  the  volume  conversion 
factor  (see  Section  IV  B) .  Normally,  the  energy  of  formation  calculated 
for  the  composition  is  used;  if  you  want  to  use  a  different  one,  it  is 
entered  as  e. 

Example: 

EXP,  V,  S 

EX1,  P,  p,  T,  t 

This  instruction  is  the  same  as  EXP,  except  that  it  is  for  a  gaseous 
mixture  at  pressure  p  and  temperature  t.  The  ideal  gas  law  is  used  to 
calculate  the  volume  and  then  EXP  is  called  automatically.  The  initial 
energy  cannot  be  changed. 

Example: 

EX1 ,  P,  2,  T,  298.15 

FLOor,  n 

This  instruction  sets  the  level  of  the  floor  at  n.  See  the  initial 
discussion  in  Section  IV.  Sometimes,  when  you  are  having  a  lot  of  trouble 
forcing  a  particular  computation  to  converge,  it  helps  to  raise  the  floor 
(that  is,  make  n  a  smaller  number  than  the  program  value).  Sometimes  it 
doesn't  help.  Once  a  new  floor  has  been  set,  it  remains  at  that  value 
until  either  a  new  COM  or  a  new  FLO  instruction  is  entered. 

FORmula,  name,  enth,  elel,  nl,  ele2,  n2,  ele3,  n3,  ... 

This  instruction  defines  'name'  as  the  name  of  a  reactant  species; 
its  enthalpy  is  'enth'  and  its  formula  contains  nl  atoms  of  element  'elel', 
n2  atoms  of  'ele2',  etc.  The  symbols  used  for  the  elements  must  be  exactly 
those  that  appeared  on  ELEment  cards  when  the  library  was  formed.  The 
enthalpy  of  formation  is  entered  in  calories  per  gram-mole,  not  calories 
per  100  grams,  and  not  kilocalories  per  mole. 

In  keeping  with  the  TIGER  and  BLAKE  conventions  for  putting  numbers 
in  instructions,  the  numbers  nl,  n2,  ...,  may  be  typed  in  any  way.  The 
program  uses  them  as  real  numbers,  so  decimal  fractions  are  accepted  and 
used  correctly.  When  the  program  goes  to  print  the  formula,  however. 
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these  numbers  are  converted  to  integer  form  by  truncation,  so  that  decimal 
fractions  will  not  be  printed  properly,  even  though  the  program  used  them 
correctly.  This  causes  trouble  with  polymers,  since  their  formulas  rarely 
contain  integer  numbers  of  atoms.  It  is  customary  for  polymer  compositions 
to  be  specified  in  units  like  gram-atoms  per  hundred  grams,  and  with  a 
corresponding  enthalpy,  but  BLAKE  does  not  accept  such  input  data.  Such 
data  are  easily  converted  to  molecular  formulas,  but  the  numbers  of  atoms 
generally  come  out  to  be  decimal  fractions.  This  creates  a  small  but 
annoying  problem,  since  one  usually  prefers  to  have  the  correct  formulas 
printed  for  the  input.  The  simplest  solution  is  to  multiply  the  formulas 
and  the  associated  enthalpy  by  either  1000  or  10000. 

For  y  r  convenience,  the  formulas  of  many  common  ingredients  of 
US  propellants  are  stored  in  the  program,  and  need  not  be  entered  on  a 
FORmula  card.  The  names  and  abbreviations  of  these  substances  are  listed 
in  Section  V. 

If  a  FORmula  instruction  is  mistyped,  or  if  a  duplicate  formula  name 
is  entered,  the  instruction  is  rejected  and  an  error  message  is  printed. 

If  you  wish  to  use  a  different  enthalpy  for  a  substance  whose  formula 
is  prestorei,  you  must  enter  a  FORmula  instruction  for  the  material  with 
a  different  name;  the  same  name  cannot  be  redefined  without  changing 
the  program  and  then  recompiling  it. 

A  maximum  of  SO  additional  formulas  is  permitted;  if  you  try  to 
enter  more,  they  are  rejected,  and  an  error  message  is  printed. 

Examples: 

F0R,  DCS,  -280E3 ,  C,  18,  H,  34,  0,  4 
F0R,  PG,  -4E7 ,  C,  4000,  H,  5000,  0,  100 


FREeze 


Most  calculations  are  carried  out  with  the  assumption  of  total 
dynamic  chemical  equilibrium.  The  FREeze  instruction  'freezes'  all 
chemistry,  so  that  all  following  calculations  are  done  on  a  mixture  of 
non-reacting  gases.  See  also  MELt. 

GEO,  name 


Change  the  Gaseous  Equation  Of  state 
Equation  of  State 
ideal 

truncated  virial 
B-K-W 
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to  one  of  the  following: 
Keyword 
IDEAL 

VIRIAL  or  LJ 
BKW 


The  default  equation  is  the  truncated  virial  equation.  If  ’name'  is 
mistyped,  the  program  defaults  to  the  truncated  virial. 

Example: 

GEO,  VIRIAL 

GRId,  CENTER,  P,  np,  ip,  np,  V,  vl,  nv,  v2 
or 

GRId,  CORNER,  P,  pi,  ip,  p2,  V,  vl,  iv,  v2 

This  instruction  computes  chemical  equilibrium  at  a  grid  of  points 
centered  about,  or  cornered  at,  a  specific  point.  The  computations  may 
be  done  for  various  combinations  of  thermodynamic  variables: 

pressure,  volume  (P,  V); 
pressure,  entropy  (P,  S) ; 
volume,  temperature  (V,  T); 
volume,  energy  (V,  E) ;  and 
volume,  entropy  (V,  S) . 

In  each  of  these  cases,  the  parenthesized  letters  are  the  key  letters 
that  must  be  used  in  the  instruction.  In  all  cases  the  order  must  be 
that  given  here. 

The  units  of  the  extensive  variables  (enthalpy,  energy,  volume, 
and  entropy)  are  specific;  that  is,  for  one  gram  of  substance. 

'pi'  is  the  initial  value  of  the  variable  (P,  in  this  case);  'p2' 
is  its  final  value.  For  CORNER,  'ip'  and  'iv'  are  the  sizes  of  the  in¬ 
crements  in  those  variables.  For  CENTER,  'np'  and  'nv'  are  the  number 
of  increments  to  be  inserted  between  initial  and  final  values.  The  in¬ 
crements  'ip'  and  'iv'  may  be  negative  if  desired. 

If  instead  of  a  number  'pi',  'p2',  etc.,  a  comma  (with  at  most 
intervening  blanks)  is  used,  the  program  will  automatically  use  the 
existing  value  of  that  variable.  See  Test  Case  Nine  for  an  example. 


Examples: 


GRID,  C0RNER,  P,  1,  0.1, 

10, 

V, 

20, 

.5, 

40 

GRID,  CENTER,  P,  1,  4, 

10, 

V, 

20, 

19, 

40 
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GUN,  s,  i,  f 


This  instruction  computes  the  thermodynamic  properties  of  the  chamber 
gas  produced  by  propellant  burning  to  equilibrium  at  constant  volume 
at  loading  densities  between  s  and  f,  at  increments  in  loading  density 
of  i.  A  maximum  of  16  different  loading  densities  is  permitted.  If  you 
try  to  use  more  than  16,  only  the  first  16  are  computed,  and  the  program 
continues  to  the  next  instruction. 

Example: 

GUN,  0.05,  .05,  0.40 
ISOline,  T,  tl,  P,  pi,  np,  p2,  (LOG) 

This  instruction  computes  thermodynamic  equilibrium  for  a  set  of 
points,  one  of  which  is  held  constant.  In  addition  to  temperature  and 
pressure  (as  shown) ,  the  same  sets  of  variables  permitted  in  a  GRId  in¬ 
struction  are  accepted. 

The  variable  to  be  held  constant  is  named  first,  'pi'  is  the  initial 
value  of  the  variable  being  changed,  'np'  is  the  number  of  increments, 
and  ' p2 '  is  the  final  value.  If  the  optional  parameter  LOG  appears, 
logarithmic  spacing  is  used.  Logarithmic  spacing  is  automatically  used 
when  the  variables  are  (S,  V). 


Example: 

ISO,  T,  1000,  P,  10,  4,  20 
LJP,  namel,  T,  vl,  S,  v2,  name2,  T,  v3,  S,  v4,  ... 

This  instruction  changes  the  values  of  the  Lennard-Jones  parameters 
for  substances  'name  1',  'name  2*,  etc.  Of  course,  'name  1'  must  be 
a  legal  name  (one  that  appears  in  the  library  as  a  CONstituent) .  The 
program  checks  that  both  parameters  have  been  entered,  even  if  you  only 
wish  to  change  one  of  them.  Changing  only  one  is  simplified  by  the 
fact  that  if  two  successive  commas  with  at  most  blanks  between  them  appear 
in  place  of  'vl'  (or  'v2',  etc.)  the  program  will  use  the  existing  value. 

The  order  of  T  and  S  is  immaterial,  but  'vl'  and  'v2'  must  correspond 
to  the  key  letters  used.  The  key  letter  T  may  be  replaced  by  E;  and  S, 
by  D. 


If  'namel'  is  mistyped  or  is  an  illegal  name,  no  action  is  taken 
but  a  diagnostic  message  is  printed.  If  two'  or  more  names  appear  on  a 
card  and  one  of  them  is  wrong,  the  program  will  nevertheless  attempt  to 
process  the  remaining  names. 
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If  no  COMposition  card  has  been  entered,  the  program  will  reject 
an  LJP  instruction.  A  suitable  message  is  printed  to  add  insult. 

Examples: 

LJP,  H20,  T,  600,  D  , 

LJP,  NH3 ,  D,  4,  T,  350 

MELt 

This  instruction  undoes  the  effect  of  a  FREeze  instruction.  It,  too, 
is  optional. 

ORC,  namel,  name2,  name3,  ... 

This  instruction  is  almost  identical  with  the  ORDer  instruction, 
except  for  one  crucial  point.  ORDer  applies  only  to  gaseous  constituents, 
while  ORC  applies  only  to  condensed  ones.  It  is  not  often  needed,  but  in 
some  cases  (e.g.,  black  powder),  it  is  quite  useful. 

Example: 

ORC,  K2S04$ 

ORDer,  namel,  name2,  name3,  ... 

This  instruction  makes  the  program  chose  ’namel',  'name2',  etc., 
as  the  components  for  the  first  attempt  at  finding  the  equilibrium  com¬ 
position.  If  these  components  are  not  all  linearly  independent,  one  or 
more  of  them  will  be  omitted.  The  names  of  the  constituents  must  be 
exactly  the  same  as  the  names  used  in  the  CONstituent  instruct’ on  when 
the  binary  library  was  formed.  Also,  none  of  these  substances  may  have 
been  previously  rejected  by  a  REJect  instruction.  If  either  of  these 
conditions  is  violated,  the  entire  case  is  omitted  and  the  program  con¬ 
tinues. 

Example: 

ORD,  N2,  H20,  CO, 

POInt,  P,  p,  T,  t 

This  instruction  computes  chemical  equilbrium  at  the  point  with  the 
specified  pressure  ’p’  and  the  specified  temperature  • t ' .  The  order  in 
which  they  appear  is  not  arbitrary;  it  must  be  (P,T),  (P,V),  (P,H),  (P,S), 
(V,E),  (V,T),  or  (V,S).  If  the  wrong  order  is  used,  an  error  message  is 
printed  and  the  program  proceeds  to  the  next  case. 
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Example: 

POINT,  P,  10,  T,  2509 

PRLimit ,  CON,  vl,  FIT,  v2,  ERR,  v3,  PAG,  v4 

This  instruction  controls  the  amount  of  output  printed  by  BLAKE. 

The  order  of  the  keywords  is  immaterial  but  they  must  be  typed  or  punched 
just  as  shown.  (For  example,  do  net  use  'PAGE'  for  'PAG'.)  The  numbers 
'vl',  'v2',  ...,  are  either  0,  1,  and  have  the  following  significance 

CON,  2:  Print  the  list  of  constants  fo'f'-the  fitting  of  the  thermo¬ 
dynamic  data  of  each  constituent.  This  list  is  printed  once  after  a  COM- 
position  instruction  is  entered  provided  PRL,  CON,  2  has  been  specified. 

The  list  is  not  repeated  until  another  PRL,  CON,  2  instruction  is  g^ven. 
Instead  the  program  reverts  to  PRL,  CON,  1. 

CON,  1:  This  is  the  default  setting.  The  list  of  constants  is  not 
printed  but  the  composition  of  the  chamber  gas  is  printed  after  each  GUN 
instruction.  PRL,  CON,  2  reverts  to  this  setting. 

CON,  0:  Tne  printing  of  both  the  constants  and  the  composition  of 
the  chamber  gas  is  suppressed.  Once  set,  this  instruction  remains  in 
effect  until  changed. 

FIT,  i:  Fit  the  four  quantities,  chamber  pressure  and  temperature, 
and  covolume  and  gamma  of  the  chamber  gas,  to  a  cubic  polynomial  in  the 
loading  density.  The  fitting  is  done  pi'ovided  at  least  four  loading  den¬ 
sities  are  called  for  in  the  GUN  instruction. 

FIT,  0:  Suopress  The  fitting,  no  matter  how  man}'  loading  densities 
are  called  for.  Once  entered,  this  instruction  remains  in  effect  until 
changed.  This  is  the  default  setting. 

PAG,  1:  Put  each  distinct  part  of  the  output  from  a  GUN  instruction 
on  a  separate  page.  These  parts  are  the  composition  (including  the  list 
of  thermodynamic  constants  if  called  for),  the  thermodynamic  functions 
for  the  chamber  gas,  and  the  summary.  This  is  the  default  setting.  Each 
output  page  produced  has  its  own  title  (if  one  was  entered) ,  the  date 
(either  from  the  system  or  from  a  DATe  instruction),  and  the  version  de- 
si  gnator . 

PRL,  0:  Do  not  paginate  the  output  until  a  page  is  filled,  or  until 
a  new  COMposition  instruction  is  entered. 

ERR,  1:  Print  each  occurrence  of  error  messages  28  and  30  as  they 
occur. 
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ERR,  0:  Suppress  the  priming  of  these  error  messages.  A  summary 
oi  the  total  number  of  occurrences  of  each  error  is  printed  at  the  end 
of  the  case.  This  is  the  default  setting. 

.'ot  all  of  the  parameters  on  a  PRL  instruction  need  be  specified 
at  once;  you  need  only  include  the  one(s)  you  want  changed. 

Examples: 

PRL,  CON,  0 
PRL,  CON,  0,  PAG.  0 
PRL,  CON,  2,  FIT,  1 

QEXplosion,  R,  r,  T,  t 

This  instruction  computes  the  heat  of  explosion  of  a  composition  at 
loading  density  *r*  at  five  freeze-out  temperatures  centered  around  tem¬ 
perature  '  t ' . 

Tne  composition  must  have  been  entered  previously.  The  order  of 
'r'  and  ' t '  is  unimportant,  but  the  values  entered  for  'r*  and  't '  must 
correspond  to  their  key  letters. 

The  five  temperatures  are  t  -  200,  t  -  100,  t,  t  +  100,  and  t  +  200. 

The  algorithm  used  is  that  of  Corner* .  The  program  executes  a  POInt 
computation  at  the  volume  l/i  and  the  chosen  temperatures. 

The  working  of  the  algorithm  is  critically  der  ndent  upon  the  contents 
of  two  DATA  statements  in  the  subroutine  that  cont:''  the  names  and  enthal¬ 
pies  of  formation  of  stable  final  products  at  298  K.  This  list  has  beer, 
selected  for  compounds  formed  by  the  elements  carbon,  hydrogen,  nitrogen, 
oxygen,  boron,  potassium,  and  sulfur.  If  additional  elements  are  present 
in  a  given  formulation,  appropriate  changes  must  be  made  to  the  DATA 
statements.  Version  295.0  allows  room  for  seven  additional  compounds 
without  ether  programming  changes  being  required. 

Example: 

QEXP,  R,  0.2,  T,  1100 

quit 

This  instruction  should  be  the  final  one  in  every  BLAKE  run.  It 
causes  the  program  to  execute  a  normal  FORTRAN  STOP.  If  it  is  omitted, 
the  program  will  likely  stop  because  it  will  next  read  an  end-of-file. 

Using  the  QUIt  instruction  is  more  elegant. 
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RECall,  namel 

Example: 

RECALL,  NUM1 

This  instruction  recalls  the  final  conditions  associated  with  the 
state  designated  by  'namel'.  If  the  point  has  not  been  saved  (see  SAVe) , 
an  error  message  is  printed,  and  succeeding  calculations  are  either  skip¬ 
ped  or  done  erroneously,  until  a  new  COMposition  instruction  is  entered. 

REJect,  namel,  name2,  name3,  ... 

This  instruction  eliminates  from  further  consideration  the  constit¬ 
uents  whose  legal  names  are  'namel',  'name2',  etc.  If  an  illegal  or 
undefined  name  is  entered,  it  is  ignored  without  harm.  See  also  DEStroy 
and  RETain.  The  maximum  number  of  constituents  that  mav  be  REJected  is 
99. 

Example: 

REJ,  C(s) ,  C2H,  C,N,  C,  CH 

RETain 

This  instruction  causes  the  program  to  retain  a  list  of  REJects  from 
one  COMposition  to  another.  If  it  is  not  used,  an  existing  REJect  list 
is  automatically  cancelled  when  a  new  COMposition  instruction  is  entered. 
The  effect  of  RETain  is  negated  by  DEStroy. 

SAVe,  name 

BLAKE  and  TIGER  start  new  thermodynamic  calculations  from  the  results 
of  a  previous  one,  except  when  initialized  by  a  COMposition  instruction. 
Then  the  program  starts  from  a  pre-set  guess  that  may  be  far  from  the 
final  composition.  The  convergence  of  a  calculation  can  often  be  made 
more  efficient  by  having  the  program  start  from  the  results  of  an  earlier 
calculation.  This  intent  must  be  made  known  in  advance  by  use  using  the 
SAVe  instruction.  The  conditions  (including  the  final  concentrations) 
just  computed  are  saved  in  memory,  and  the  name  'name'  is  associated  with 
them.  Up  to  8  points  may  be  saved  in  this  way,  and  later  recalled  (see 
RECall).  If  you  try  to  save  more  than  8,  an  error  message  is  printed, 
and  the  point  is  not  saved,  A  name  may  be  reused,  however,  the  new  point 
will  then  overwrite  the  old  one. 


Example: 
SAV,  NUM1 


WIW WWJPWJIWH! 


STArt  the  Library 


This  instruction  is  used  to  form  the  binary  library.  See  Section 


TITle,  (message) 

This  instruction  causes  a  title  to  be  written  at  the  top  of  each  page 
of  output.  All  of  the  characters  in  ‘message'  are  centered  and  printed. 

Example: 


TIT,  THIS  IS  A  TITLE 


TIMe,  t 


This  instruction  is  system-dependent  (see  Section  VA1).  If  you  have 
implemented  it  on  your  system,  then  it  is  used  to  change  the  default  setting 
of  the  maximum  time  allowed  for  a  single  calculation  by  BLAKE.  It  is  use¬ 
ful  in  the  case  that  you  want  to  run  a  number  of  different  cases  that  you 
anticipate  will  require  a  maximum  of  (say)  five  minutes.  It  is  unlikely 
that  you  will  want  all  five  minutes  to  be  devoted  to  one  case.  By  setting 
the  TIMe  instruction,  you  can  control  approximately  the  amount  of  time 
the  program  devotes  to  one  case  before  going  on  to  the  next.  This  time 
limit  is  different  from,  and  independent  of,  the  time  limit  that  is  usually 
set  at  the  beginning  of  a  run  by  your  job  card  (or  other  job  control  entry). 

This  instruction  is  optional.  If  you  have  implemented  the  CLOCK 
instructions,  you  can  then  set  a  default  time  limit.  For  the  UNIVAC  1108, 
2.5  minutes  was  a  reasonable  value;  for  the  CDC  7600,  0.4  minutes  is  better. 

The  actual  time  that  elapses  before  time  is  called  by  this  instruction 
is  usually  greater  than  the  value  set,  because  the  program  only  checks  the 
clock  at  certain  cycle  limits. 

Example: 

TIM,  0.1 

TRAn sport  properties 

This  instruction  is  inactive.  Someday  it  will  cause  the  program  to 
compute  approximate  values  for  viscosity  and  thermal  conductivity.  Some 
day  gold  will  fall  from  the  sky. 


UN Its,  name 

This  instruction  causes  the  program  to  include  or  exclude  English 
(or  conventional  engineering)  units  for  the  summary  output  in  a  GUN 
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calculation.  The  summary  output  is  always  printed  in  SI  units.  This 
cannot  be  suppressed  except  by  changing  the  program.  If  in  addition, 
you  want  English  units  printed,  this  instruction  will  do  that  when  'name' 
is  ENG.  If  'name'  is  'SI',  the  English  unit  print  is  suppressed. 

Example: 

UNI,  ENG 

E.  Error  Messages 

There  are  two  score  plus  four  error  messages  contained  in  BLAKE. 

Only  a  few  of  them  (e.g.,  types  9  and  24)  indicate  truly  serious  problems 
and  terminate  a  case.  Most  of  the  others  call  your  attention  to  situations 
that  may  need  corrective  action,  but  the  program  continues  with  the  same 
case.  A  few  indicate  that  the  program  is  having  problems  reaching  equil¬ 
ibrium.  Where  appropriate  the  program  prints  the  instruction  that  is  the 
source  of  a  problem. 

All  of  the  error  messages  incorporated  in  the  code  are  listed  here. 
(Some  of  the  error  messages  of  TIGER  were  removed,  hence  the  gaps  in  the 
numbering.)  The  exact  message  is  listed  IN  CAPITAL  LETTERS;  then  the 
probable  cause  is  discussed  and  remedial  action  suggested. 

Error  Type  1:  INPUT  CARD  NOT  ACCEPTABLE. 

The  instruction  entered  is  not  one  of  the  37  recognized  by  BLAKE. 

The  instruction  is  skipped  and  the  program  attempts  to  interpret  and  follow 
the  next  instruction.  This  message  usually  arises  from  faulty  typing  or 
punching. 

Error  Type  2;  TOO  MANY  CONSTITUENTS  REJECTED. 

Since  up  to  99  constituents  may  be  rejected,  this  message  is  very 
improbable.  It  could  occur  towards  the  end  of  a  very  long  run  in  which 
you  forgot  to  DEStory  several  previous  REJect  lists;  even  this  is  im¬ 
probable. 

Error  Type  3:  TOO  MANY  CONSTITUENTS  ORDERED. 

Only  25  constitutents  may  be  ordered.  The  program  continues  but 
only  the  first  25  constituents  are  ordered.  There  is  never  any  need  to 
order  as  many  as  10  constituents;  one  per  element  is  sufficient. 

Error  Type  4:  TOO  MANY  CONSTITUENTS  CHOSEN. 

The  limit  is  25.  The  program  continues  but  only  the  first  25  con¬ 
stituents  are  chosen.  Again,  there  is  little  reason  for  trying  to  choose 
26  or  more  constituents.  As  mentioned  earlier,  this  order  is  better  not 
used. 
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Error  Type  6:  FINDING  RH0O  EXCESSIVE  ITERATIONS. 

This  error  message  arises  in  Subroutine  FINDRH;  PO  and  P  are  printed 
at  the  same  time.  When  P  =  PO,  a  value  for  RHO  has  been  found,  so  the 
difference  between  the  two  is  a  measure  of  the  error  at  that  time.  The 
program  uses  the  last  value  of  RHO  and  continues.  If  it  converges  later, 
this  error  can  be  safely  ignored. 

Error  Type  7:  FINDING  TO  EXCESSIVE  ITERATIONS. 

This  message  arises  in  Subroutine  F INDTP;  PO  and  P  are  printed  at  the 
same  time.  When  these  two  are  equal,  a  value  for  TO  has  been  found  so 
their  difference  indicates  the  magnitude  of  the  error.  The  program  uses 
the  last  value  of  TO  and  continues.  If  it  later  converges,  this  message 
can  be  ignored. 

Error  Type  8:  COMPONENTS  LINEARLY  DEPENDENT. 

This  message  arises  in  Subroutine  SELECT.  It  usually  means  that 
you  have  CHOsen  constituents  unwisely  or  have  been  unlucky  with  an  ORDer 
instruction.  REJect  a  condensed  constituent  and  try  it  again. 

Error  Type  9:  NO  ACCEPTABLE  COMPONENTS. 

This  is  one  of  the  most  frustrating  messages  of  all.  It  arises  in 
Subroutine  SELECT  and  indicates  that  the  program  cannot  find  a  set  of 
linearly  independent  components  after  trying  many  sets  above  the  floor. 

The  case  is  rejected  and  the  program  continues.  Unfortunately  it  does 
not  give  a  list  of  the  constituents  it  has  considered  nor  the  position 
of  the  floor.  There  is  no  certain  procedure  to  follow  at  this  point. 

I  usually  add  another  REJect  list,  play  with  an  ORDer  or  ORC  instruction, 
or  as  the  last  resort,  raise  the  floor. 

Error  Type  11;  EXPLOSION  STATE  CARD  NOT  ACCEPTABLE. 

This  message  arises  in  Subroutine  EXPLOS.  It  indicates  that  either 
an  EXPlosion  or  E/X  instruction  has  been  incorrectly  punched.  The 
offending  instruction  is  printed  to  facilitate  the  detection  of  the  error. 
The  case  is  rejected  and  the  program  continues. 

Error  Type  12:  INPUT  TO  POINT  NOT  ACCEPTABLE. 

This  message  arises  in  Subroutine  POINNT.  It  indicates  that  an  error 
has  been  made  on  either  a  POInt,  ISOline,  or  GRId  instruction.  Generally 
this  means  that  the  variables  were  ordered  incorrectly;  e.g.,  V,  P  instead 
of  P,  V.  The  instruction  is  rejected  and  the.  program  continues. 
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Error  Type  13:  THIS  POINT  NOT  SAVED. 

This  message  is  printed  in  Subroutine  RECALL  and  is  self-explanatory. 
Usually  it  arises  from  a  punching  or  typing  error-  It  can  also  occur  when 
a  previous  case  was  rejected  or  failed  to  converge,  so  the  point  that  was 
to  have  been  saved  was  not  computed.  The  program  continues  with  the 
succeeding  instruction.  The  remedy  is  obvious. 

Error  Type  14:  TOO  MANY  POINTS  SAVED. 

This  message  arises  in  Subroutine  SAVE  and  is  self-explanatory.  The 
limit  is  eight.  The  point  is  not  saved,  the  program  continues,  and  will 
probably  soon  come  to  grief  when  you  try  to  RECall  the  unsaved  point. 

The  remedy  is  obvious. 

Error  Type  16:  LIBRARY  CARD  NOT  ACCEPTABLE. 

This  message  arises  in  Subroutine  LIBRARY.  It  indicates  that  a 
library  instruction  has  failed  to  pass  all  of  the  checks  built  into  the 
program.  The  instruction  is  printed;  it  should  be  compared  with  the 
prescribed  format  for  library  instructions.  Punching  or  typing  errors 
are  the  most  common  causes.  The  program  ignores  the  bad  instruction  and 
continues,  but  you  probably  will  not  want  to  use  the  library  thus  formed. 
The  binary  library  cannot  be  corrected;  it  must  be  formed  all  at  once. 

Fix  the  card  and  start  again. 

Error  Type  17:  COMPOSITION  CARD  NOT  ACCEPTABLE. 

This  message  arises  in  Subroutine  COMPOS •  Its  most  common  cause  is 
that  a  formula  card  was  either  omitted  entirely,  or  that  it  was  previously 
rejected  by  the  program.  Punching  or  typing  errors  are  the  next  most 
common  causes.  The  offending  card  or  instruction  is  printed.  The  pro¬ 
gram  ignores  the  instruction  and  tries  to  continue  since  later  COMposi-tion 
instructions  may  be  acceptable.  The  remedy  it  obvious. 

Error  Type  IS:  FORMULA  ELEMENTS  NOT  IN  LIBRARY. 

This  message  arises  in  Subroutine  FORMUgA,  and  is  self-explanatory. 
The  program  ignores  the  instruction  and  continues.  If  not  due  to  a  mere 
typing  blunder,  this  error  is  difficult  to  £jx,  since  it  involves  adding 
another  element  to  the  library,  a  time-consuming  job. 

Error  Type  20:  NOT  ALL  ORDERED  CONSTITUENTS  IN  LIBRARY. 

This  message  arises  in  Subroutine  COMPOS*  It  is  caused  by  the  list¬ 
ing  on  an  ORDer  or  ORC  instruction  of  a  constituent  that  the  program  can¬ 
not  find  in  the  binary  library.  Typing  or  punching  errors  are  the  first 
thing  to  look  for,  but  an  equally  likely  ertut  is  that  the  substance  in 
question  had  been  previously  listed  in  a  REJect  instruction.  The  program 
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ignores  the  instruction,  rejects  the  case,  and  continues.  The  remedy  is 
obvious  although  it  is  not  always  easy  to  decide  which  constituent  in  an 
ORDer  or  ORC  instruction  was  at  fault. 

Error  Types  21,  22  and  23:  FREEZE  MATRIX  SINGULAR. 

FREEZE  FORCE  ITERATION  EXCEEDED. 

FREEZE  ITERATION  EXCEEDED. 

All  of  these  messages  arise  in  Subroutine  FREEZE.  All  three  errors 
are  fatal  to  the  computation  being  attempted.  The  program  skips  it  and 
looks  for  a  new  COMposition  instruction.  The  usual  causes  are  either  the 
frozen  composition  requested  cannot  exist  at  the  point  requested,  or  the 
point  in  question  is  too  far  from  the  last  computed  point. 

Error  Type  24:  THERMO  ITERATION  EXCEEDED. 

This  message  arises  in  Subroutine  THERMO  and  is  fatal.  The  program 
continues  and  looks  for  a  new  COMposition  instruction.  The  most  common 
cause  is  that  it  was  not  working  with  a  good  set  of  components.  (Read 
Section  VIIA  over).  The  only  cure  is  to  use  ORDer  or  ORC  instructions 
to  rearrange  the  order  in  which  components  are  selected,  or  to  change 
the  REJect  instructions,  or  both.  It  is  worth  trying  extreme  measures 
to  make  a  calculation  converge  somehow;  once  it  has,  that  point  can  be 
SAVed  and  used  to  restart  a  new  calculation  with  fewer  restrictions. 

Error  Type  25:  MELTING  TEMPERATURE  EXCESSIVE  ITERATIONS. 

This  message  arises  in  Subroutine  ECOMPO.  It  is  not  a  fatal  error 
but  indicates  that  the  program  is  having  some  difficulty  converging. 

Keep  in  mind  that  the  version  of  TIGER  on  which  BLAKE  is  based  often 
encountered  insurmountable  problems  when  it  tried  to  deal  with  a  con¬ 
stituent  in  the  neighborhood  of  its  melting  point.  The  new  TIGER  handles 
such  computations  with  ease.  If  the  computation  you  are  trying  eventually 
converges,  fine;  otherwise  it  is  beyond  BLAKE's  capabilities. 

Error  Type  26:  ECOMPO  FORCE  EXCESSIVE  ITERATIONS. 

This  is  a  non-fatal  error  arising  in  Subroutine  ECOMPO.  It  indicates 
that  the  program  is  having  trouble  converging.  If  it  eventually  does 
converge  there  is  no  problem;  if  not,  you  will  get  a  different  message 
(usually  No.  24). 

Error  Type  27:  SOLVING  FOR  ETAS  MATRIX  SINGULAR. 

See  No.  26. 

Error  Type  28:  SOLVING  FOR  ETAS  EXCESSIVE  ITERATIONS. 

See  No.  26. 
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Error  Type  29:  SOLVING  FOR  N(I)S  MATRIX  SINGULAR. 

See  No.  26. 

Error  Type  30:  SOLVING  FOR  N(I)S  EXCESSIVE  ITERATIONS. 

See  No.  26.  Messages  28  and  30  are  often  produced  during  the  course 
of  a  straightforward  computation  that  converges  in  a  few  more  iterations. 
For  this  reason  BLAKE  generally  suppresses  their  printing.  If  you  want 
to  see  them,  use  the  instruction,  PRL,  ERR,  1.  Note  that  even  if  they 
are  suppressed,  a  summary  of  the  number  of  times  each  was  encountered 
is  still  printed. 

It  is  of  course  possible  to  make  the  appearance  of  all  these  error 
messages  less  likely  by  raising  the  limits  of  the  iteration  counters  set 
in  the  main  program.  I  have  played  this  game  occasionally,  but  while 
it  always  used  up  more  computer  time,  it  rarely  helped  with  a  difficult 
computation. 

Error  Type  31:  TOO  MANY  GASEOUS  CONSTITUENTS. 

This  message  arises  in  Subroutine  COMPOS,  and  indicates  that  more 
than  29  gaseous  constituents  have  been  selected  from  the  binary  library. 

It  does  not  tell  you  which  29  they  were  but  you  can  readily  work  this 
out  for  yourself  by  reference  to  a  library  print.  The  program  skips  this 
case  and  continues  looking  for  a  new  COMposition  instruction.  The  only 
cure  is  to  add  more  REJects. 

Error  Type  32:  TOO  MANY  CONDENSED  CONSTITUENTS. 

The  limit  here  is  10,  but  experience  shows  that  the  program  will 
not  converge  when  more  than  four  condensed  constituents  have  been  selected. 
With  the  library  furnished  with  BLAKE  this  error  can  arise  only  if  you 
try  compositions  with  at  least  nine  elements. 

Error  Types  33  -  37:  TOO  MANY  (  _ _  )  CARDS. 

All  of  these  errors  arise  in  Subroutine  LIBRARY  and  are  self-explan¬ 
atory.  The  ellipsis  can  be  filled  with  ELEMENTS,  CONSTITUENTS,  STG,  STR, 
or  STC.  The  limits  presently  set  in  the  program  are: 

ELEments:  40 

CONstituents:  180 

STG:  99 

STR:  S40 

STC:  180 
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Any  of  these  errors  can  be  readily  avoided  by  redimensioning  the  ap¬ 
propriate  array  in  Subroutine  LI  BRAY.  If  you  do,  remember  to  change  the 
error  messages  accordingly.  (This  is  not  absolutely  necessary,  of  course; 
it  is  likely  that  the  error  message  will  be  encountered  by  someone  else 
who  will  have  no  idea  that  it  is  meaningless.  As  booby  traps  go,  this  is 
a  relatively  harmless  one.) 

The  last  six  error  messages  unnumbered: 

TWO  OR  MORE  (...)  HAVE  THE  SAME  NAME 

The  ellipsis  can  be  either  ELEMENTS  or  CONSTITUENTS.  Both  messages 
arise  in  Subroutine  LIBRAY  and  are  self-explanatory.  If  you  had  to  limit 
the  length  of  a  name  to  four  characters  (see  Section  VI  A3),  you  will 
surely  see  this  message  unless  you  have  scrupulously  examined  the  library 
before  attempting  to  form  the  binary  library  file  for  the  first  time. 

The  remedy  is  obvious. 

LJP  CARD  IS  FAULTY  OR  PREMATURE. 

This  message  arises  in  Subroutine  LJCH  and  is  self-explanatory.  The 
program  skips  the  instruction  and  proceeds. 

---DOES  NOT  APPEAR  IN  THE  LIST  OF  GASEOUS  CONSTITUENTS 

The  dash  is  filled  by  the  name  of  a  putative  constituent  from  a  LJP 
instruction  (see  Section  VII  C2) .  The  meaning  is  self-explanatory.  The 
program  skips  the  instruction  and  proceeds. 

NO  COMP  CARD,  NO  CALCULATION 

This  message  arises  in  Subroutine  QEXP  and  is  self-explanatory.  The 
program  skips  the  instruction  and  proceeds. 

CALCULATION  REJECTED 

This  message  arises  in  Subroutine  REJECT.  When  either  fatal  error 
9  or  24  is  encountered,  the  program  prints  the  appropriate  error  message 
and  then  branches  to  Subroutine  REJECT  where  various  counters  are  reset 
and  this  message  is  printed.  The  program  then  proceeds,  but  almost  all 
further  instructions  are  skipped  until  a  new  COMposition  instruction  is 
found.  Only  TITie  and  QUIt  are  recognized. 

It  is  hard  to  give  good  advice  at  this  point;  however,  see  the  discus¬ 
sion  for  Error  Types  9  and  24. 

UNITS  CARD  IS  FAULTY.  NOTHING  DOME. 
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This  message  arises  in  the  main  program  and  is  self-explanatory. 

Only  ENG  or  SI  are  recognized  keywords  in  a  UNIts  instruction.  The  pro¬ 
gram  proceeds. 

GEO  CARD  IS  FAULTY. 

This  message  arises  in  the  main  program.  Only  IDEAL,  L-J,  BKW,  and 
VIRIAL  are  recognized  keywords  on  a  GEO  card.  If  a  faulty  card  is  entered, 
the  program  reverts  to  the  truncated  virial  equation  of  state,  and  then 
proceeds. 

ILLEGAL  INPUT  FOR  CJPREP.  OMITTED. 

This  message  arises  in  Subroutine  CJPREP  and  indicates  that  an  error 
was  made  on  an  EX1  instruction.  The  instruction  is  skipped  and  the  program 
proceeds  with  the  next  instruction. 

TIME  HAS  BEEN  CALLED.  RETURN  TO  START.  DO  NOT  PASS  'GO*.  DO  NOT 

COLLECT  TWO  HUNDRED  DOLLARS. 

As  mentioned  in  Section  VI  Al,  the  program  contains  the  capability 
of  timing  a  calculation  independent  of  the  total  time  allowed  for  the 
entire  run,  provided  you  have  implemented  the  call  to  the  system  clock. 

This  message  is  printed  if  the  allotted  time  is  exceeded.  There  are  few 
calculations  that  will  converge  in  a  longer  time  if  the  default  time  is 
appropriate  to  your  computer.  (The  BRL's  experience  has  been  that  0.4 
minute  is  a  suitable  limit  for  the  CDC  7600 ,r  and  2  minutes  for  a  UNIVAC 
1108.  A  few  extra  lines  of  programming  will  readily  enable  you  to  print 
out  the  elapsed  time  per  run  if  you  have  previously  implemented  the  call 
to  the  system  clock.) 

Rather  than  setting  a  longer  time  limit  (using  the  instruction  TIMe, 
t) ,  you  should  fiddle  with  the  ORDer  and  REJect  instructions  to  try  to 
make  the  calculation  converge  more  rapidly. 

DUPLICATE  FORMULA  CARD  --  IGNORED. 

This  message  arises  in  Subroutine  FORMLA;  it  is  followed  by  a  listing 
of  the  duplicate  card  in  question.  It  means  that  you^have  entered  a  FORmula 
card  using  a  reserved  name  on  the  prestored  list.  The  program  checks 
only  on  the  name,  not  on  the  enthalpy.  If  you  want  to  use  a  different 
enthalpy  for  a  pre-stored  formula,  you  must  insert  a  new  FORmula  card 
with  a  different  name. 

BAD  FORMULA  CARD  --  EITHER  TOO  MANY  OF  THEM  OR  THE  CARD  IS  MISPUNCHED. 

This  message  arises  in  Subroutine  FORMLA  and  is  self-explanatory. 

In  checking  to  see  why  the  program  has  declared  a  card  to  be  mispunched, 
pay  particular  attention  to  the  difference  between  letter  0  and  numeral 
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0.  In  any  case  the  offending  card  is  skipped  and  the  program  proceeds. 

Soon  afterwards,  you  will  likely  get  Error  Type  17.  Such  is  life. 

ERROR.  GUN  CARD  HAS  MIS-FIRED. 

This  message  arises  in  Subroutine  GUNCAL  and  means  that  a  GUN  instruc¬ 
tion  is  faulty.  The  program  skips  it  and  proceeds. 

AN  'EOF'  HAS  BEEN  ENCOUNTERED.  EXECUTION  TERMINATED. 

This  message  arises  in  Subroutine  CARDRD.  On  CDC  computers  it  prevents 
the  program  from  aborting  when  it  tries  to  read  an  end-of-file,  thus  sav¬ 
ing  a  page  of  incomprehensible  printing.  You  can  always  end  the  program 
this  way  but  it  is  suggested  that  you  enter  a  QUIt  instruction  instead. 

These  44  messages  are  the  total  embedded  in  BLAKE;  any  others  must 
have  come  from  your  operating  system. 

VIII.  THE  OUTPUT  FROM  'BLAKE' 

Apppndix  C  contains  the  output  produced  by  the  BRL's  CYBER  system 
for  nine  different  BLAKE  runs.  They  are  included  here  both  for  their 
didactic  value  and  as  comparisons  for  users  to  ensure  that  their  computers 
are  producing  the  same  numbers.  In  making  such  comparisons  bear  in  mind 
that  agreement  for  the  last  digit  is  not  to  be  expected;  discrepancies 
of  up  to  a  few  tenths  of  a  percent  can  be  expected  from  rounding  errors. 

The  ECHo  feature  was  used  in  all  of  these  runs  to  show  the  instruc¬ 
tions  entered.  Special  attention  should  be  given  to  the  REJect,  ORDer, 
and  ORC  instructions  since  these  potentially  offer  the  most  difficulty 
for  users. 

The  output  starts  with  a  banner  page  identifying  the  version  number 
and  proclaiming  BLAKE'S  relation  to  the  TIGER  program.  The  date  is 
printed  next  (assuming  either  that  the  system  function  DATE  was  imple¬ 
mented,  or  that  a  DATe  instruction  was  used). 

If  either  a  CMT  or  ECHo  instruction  was  entered,  their  output  comes 
on  the  same  page;  neither  of  these  instructions  ever  generates  a  page 
eject. 

A  COMposition  instruction  always  produces  a  page  eject.  The  auto¬ 
matic  page  numbering  starts  with  this  next  page;  it  cannot  be  suppressed 
or  re-started  within  a  run.  If  a  TITle  instruction  was  entered,  the  title 
(centered)  is  printed.  Then  comes  the  standard  heading,  THE  COMPOSITION 
IS,  which  is  followed  by  a  listing  of  the  ingredients  specified  in  the 
COMposition  instruction. 
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Each  ingredient  is  followed  by  its  v/eight  and  mole  percent  in  the 
formulation,  its  enthalpy  of  formation,  and  its  formula.  In  the  test 
cases,  note  that  most  of  the  ingredients  are  those  already  stored  in  the 
program,  and  are  called  up  by  using  the  proper  abbreviation.  One  except 
tion  occurs  in  the  test  case  two  where  it  was  desired  to  use  potassium 
cryolite  as  an  ingredient  instead  of  sodium  cryolite.  The  latter  is 
stored  in  the  program  under  the  name  CRY;  potassium  cryolite  must  there-- 
fore  be  given  a  different  name  (KCRY,  in  this  case) . 

The  program  next  prints  the  heat  of  formation  of  the  composition  oh 
both  specific  (per  gram)  and  molar  bases.  (Owing  to  the  form  of  the 
formulas  used  for  nitrocellulose  and  other  polymers,  their  mole  fractions 
are  distorted;  likewise,  the  molar  heat  of  formation  of  the  formulation 
is  rarely  meaningful. 

Then  the  program  prints  the  ultimate  elemental  composition  of  the 
formulation  in  mole  percent. 

A  page  eject  normally  occurs  at  this  point  unless  it  has  been  sup.- 
pressed  by  a  PRL,  PAG,  0  instruction.  (See,  for  example,  the  output  for 
test  case  5.)  If  the  list  of  possible  products  has  been  requested  (by  a 
PRL,  CON,  2  instruction),  it  is  now  printed.  See  page  2  of  the  output 
for  test  case  one. 

After  the  standard  heading  and  title  (if  any)  are  printed,  the 
program  prints  the  thermodynamic  constants  for  the  gaseous  constituents 
of  possible  products.  Not  all  of  these  may  ultimately  appear  in  the 
final  equilibrium  composition,  but  no  others  will  be  considered.  The 
ordering  is  the  one  in  which  these  constituents  appear  in  the  binary 
library  unless  an  ORDer  or  CHOose  instruction  has  been  used. 

The  information  printed  for  each  constituent  is  as  follows.  The 
formula  of  the  constituent  (as  defined  previously  on  a  CONstituent  card 
when  the  binary  library  was  formed),  the  BKW  covolume,  the  two  Lennard- 
Jones  force  constants,  and  then  the  9  constants  that  reproduce  the  various 
thermodynamic  properties  of  the  constituent  as  a  function  of  temperature. 
(See  the  discussion  in  Section  VII  Cl,  under  STR.) 

The  floor  (see  Section  VII  A)  is  printed. 

If  any  condensed  constituents  have  been  selected,  their  thermodynamic 
data  are  printed  next.  Owing  to  the  length  of  the  names  of  some  of  the 
condensed  species,  some  arbitrary  truncation  had  to  be  imposed  on  them* 

The  correct  formulas  are  printed  in  the  listing  produced  when  the  library 
is  first  formed.  This  is  the  explanation  of  the  cryptic  message  FRMS 
IN  LIB  PRT  that  next  appears.  The  symbols  B  and  C  are  the  coefficients 
for  the  thermodynamic  data  and  the  equation  of  state  data  for  the  condensed 
species.  (See  STR  and  STC  in  Section  VII  Cl.)  As  noted  on  v.he  print,  the 
data  are  arranged  in  the  order  thermo  data  for  liquid,  then  solid,  and 
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p-V-T  data  for  liquid,  then  solid.  If  the  species  in  question  does  not 
form  a  liquid  (solid)  phase,  only  zeroes  appear  for  the  liquid  (solid) 
line. 


The  next  part  of  the  output  depends  on  the  computation  selected. 

If  the  printing  of  the  final  equilibrium  composition  has  not  been  suppres¬ 
sed,  it  comes  next.  The  order  is  now  different  from  that  in  which  the 
thermodynamic  constants  were  printed;  it  is  essentially  that  of  decreasing 
concentration.  The  concentration  unit  printed  is  moles  per  kilogram  of 
starting  materials  which  has  proved  the  most  convenient  for  users.  The 
program  prints  the  concentrations  of  all  gaseous  species  first  followed 
by  the  concentrations  of  condensed  species;  then  it  prints  the  total 
moles  of  gas.  Eventually  the  total  grams  of  condensed-phase  material 
will  be  printed  but  this  it  not  now  done  in  version  205.0. 

For  a  GUN  instruction  the  program  next  prints  a  summary,  the  outpat 
of  most  interest  to  ballisticians.  This  summary  presents  the  properties 
of  the  chamber  gas  at  equilibrium  for  each  loading  density  specified  in 
the  GUN  instruction.  These  properties  are  the  pressure,  temperature, 
impetus,  molecular  weight  of  the  gas,  covolume,  gamma,  heat  capacity  at 
constant  pressure,  the  second  and  third  virial  coefficients,  the  specific 
gas  volume,  and  also  the  entropy,  enthalpy,  internal  energy,  adiabatic 
exponent,  and  compressibility. 

Most  of  these  quantities  are  self-explanatory.  If  no  condensed 
phase  has  formed  then  the  specific  gas  volume  is  exactly  the  reciprocal 
of  the  loading  density;  if,  however,  some  condensed  material  has  formed, 
BLAKE  subtracts  its  volume  from  the  total  (geometrical)  volume  and  prints 
the  remainder.  It  is  this  remaining  volume  that  is  printed  under  the 
heading,  GAS  VOL. 

Likewise,  the  heading  MOL  WT  GAS  indicates  that  the  number  there 
printed  refers  to  the  gas  phase  only.  This  point  should  be  kept  in  mind 
when  comparing  BLAKE  output  with  that  produced  by  CEC71,  where  the  mole¬ 
cular  weight  printed  is  for  the  entire  equilbrium  mixture.  Also,  this 
correction  carries  over  to  the  IMPETUS,  which  in  BLAKE  is  computed  strict¬ 
ly  for  the  gas  phase  only. 

The  heat  capacity  and  the  gamma  printed  are  computed  with  the 
assumption  that  the  chemical  composition  of  the  mixture  is  fixed,  or 
'frozen.'  The  equilibrium  heat  capacity  is  larger,  since  it  contains 
the  frozen  term  plus  a  contribution  arising  from  the  effect  of  tempera¬ 
ture  or.  the  shifting  equilibrium.  This  same  effect  makes  the  equilibrium 
gamma  smaller.  The  difference  is  quite  noticeable  for  hot  propellants 
like  M9.  The  matter  is  discussed  in  more  detail  by  several  writers;  see, 
for  example,  a  paper  by  Fifer  and  Cohen24.  BLAKE  prints  only  the  frozen 


A.  Fifer  and  A.  Cohen ,  "High  Temperature  Thermal  Conductivity  and 
Heat  Transfer  in  Gun  Propellant  and  Primer  Combustion  Gases , "  Proceedings 
of  the  XIV  JANNAF  Combustion  Meeting. 
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quantities  since  these  are  the  ones  that  have  been  traditionally  used  in 
ballistic  performance  codes.  Whether  this  is  correct  or  not,  and  the  size 
of  the  ?rror  thus  introduced  if  it  is  not,  remain  open  questions. 

The  adiabatic  exponent  is  the  true  equilibrium  derivative  -(31nP/31nP)s. 

The  quantity  $  that  is  printed  last  is  the  ratio  cf  the  non- ideal 
pressure  to  the  ideal  pressure.  It  is  thus  a  measure  of  the  deviation  of 
the  system  from  ideality.  The  theory  on  which  the  gas  pressure  correction 
is  based  may  not  be  valid  for  $  much  greater  than  1.5. 

If  the  instruction  UNIts,  ENG,  has  been  entered,  the  program  now  repeats 
the  summary,  this  time  using  the  proscribed  units,  pounds  mass,  pounds 
force,  pounds  per  square  inch,  cubic  inches,  etc. 

If  four  or  more  loading  densities  have  been  requested  and  if  a  PRL, 

FIT,  1  instruction  has  been  entered,  the  next  part  of  the  output  con¬ 
sists  of  coefficients.  These  coefficients  are  for  cubic  polynomials  whose 
independent  variable  is  the  loading  densicy.  They  are  printed  in  order 
of  increasing  power  from  0  to  3.  The  quantities  fitted  are  the  impetus, 
chamber  pressure,  chamber  temperature,  and  covolume.  The  coefficients 
are  obtained  by  a  least-squares  process  and  are  useful  from  the  lowest 
density  up  to  the  largest  density  fitted.  They  should  not  be  used  for 
even  short  extrapolations. 

Test  cases  six,  seven,  and  eight  are  further  examples  of  the  use  of 
BLAKE  for  GUN  calculations.  Note  the  variations  produced  by  the  use  of 
different  PRL  instructions.  These  cases  also  give  interesting  examples 
of  the  application  of  the  instructions  ORDer,  REJect,  and  ORC. 

Test  case  nine  shows  a  different  calculation  that  is  sometimes  useful. 

A  typical  double-base  propellant  is  first  allowed  to  reach  equilibrium 
at  a  constant  volume  of  5  cc/g;  the  resulting  hot  gas  is  then  expanded 
adiabatically.  A  PRL  instruction  is  used  to  control  the  printing  of  the 
constituents. 


The  expansion  calculation  is  ordered  by  the  instruction 
ISOline,  S,  ,  V,  5,  7,  500 


Note  the  double  ,  ,  ;  this  instructs  the  program  to  use  the  value  of 
entropy  previously  computed,  a  number  the  user  does  not  know  in  advance. 


The  output  in  this  case  is  mainly  self-explanatory. 
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printed  here  are,  respectively,  the  quantities  — 
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SIGMA  is  the  Riemann  integral. 
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where  =  local  sound  speed,  and 


a  =  log  v 

OMEGA  is  the  integral 
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APPENDIX  A. 

COMPOSITIONS  OF  PROPELLANTS  USED  IN  COMPARISON  CALCULATIONS 


APPENDIX  A 


COMPOSITIONS  OF  PROPELLANTS  USED  IN  COMPARISON  CALCULATIONS 

II.  Vest 

Designation:  Type  I  Type  II 


Ingredient 

NC  1313 

83.25 

— 

Nitroglycerin 

10.29 

20.49 

Dibutylphthalate 

5.14 

— 

Diphenyl amine 

1.01 

0.70 

Water 

0.31 

0.41 

NC  1323 

— 

78.46 

Dinitrotoluene 

— 

0.97 

Potassium  nitrate 

— 

0.94 

Barium  nitrate 

— 

1.36 

Graphite 

— 

0.27 

74 


APPENDIX  B 

FORMULAS  OF  INCREDIENTS  PRESTORED  IN  'BLAKE 


^j?i> 


gp ->j,  — 


1 

s 

i 

I 

| 

BLAKE 

NAME 

HEAT  CF 
formation 

CAL/MOL6 

MOL. 

WEIGHT 

FORMULA 

| 

ACE  TOM 

-.59  270E+.05 

58.0804 

C  H  0 

3  6  1 

ETOH 

-.66420E+05 

46.0694 

C  H  0 

2  6  1 

i 

i 

EC 

-• 2510QE+05 

268.3598 

CHON 

17  20  1  2 

' 

s 

NDP  A 

•15400E+05 

214.2242 

CHON 

12  10  2  2 

f, 

£ 

; 

NNOP 

• 54000E+04 

259.2217 

CHON 

12  9  4  3 

: 

*' 

£ 

NG 

-.886006+05 

227.0872 

CHON 

3  5  9  3 

1 

NQ 

-.22100E+C5 

104.0686 

CHON 
14  2  4 

! 

DNT 

-• 17100E+05 

182.1360 

CHON 

7  6  4  2 

i  * 

; 

; ! 

DBP 

-• 20140E+06 

278.3496 

C  H  0 

16  22  4 

\\ 
i  \ 

\  f 

i 

ALC 

-• 66420E  +05 

46.0694 

C  H  0 

2  6  1 

[f 

KS 

-•34266E+06 

174.2536 

K  S  0 

2  14 

I 

■E  E 

'  s 

t  J 

CRY 

-.789006+06 

209.9413 

NA  AL  F 

3  16 

r  > 

\\ 

1  > 

KN 

-.117766+06 

101.1029 

K  N  0 
113 

i  ‘ 

F  ; 

t  : 

banitr 

-.23706E+06 

261.3498 

BA  N  0 
12  6 

$-  L 

^  f 
£*  « 

£  1 
|  ! 

£  ' 

DPA 

H20 

.  31|07,OE+P5 

-.603156+05 

169.2267 

18.0154 

C  N  H 

12  1  11 

H  0 

2  1 

t  \ 

i  \- 
%  =- 
a.  . 

|  :• 

‘T  t 

j  l 

?■  ?- 
& 

II 

|;  £ 

‘5 

l| 

1  i 

#  ._ 

C 

0. 

12.011P 

77 

C 

1 

77 


BLAKE 

NAME 


HEAT  OF 

FORMATION 

CAL/MOLE 


MOL. 

HEIGHT 


FORMULA 


KBF4 

-• 45100E+06 

125.9016 

K 

B  P 

1  1 

4 

PETN 

-.12880E+06 

316*1366 

C 

H  0 

N 

5  0  12 

4 

ROX 

• 14690E+05 

222.1176 

C 

H  N 

0 

3  6 

6 

6 

AMYLPH 

-.22090E+P6 

306.4036 

C 

H  0 

18  26 

4 

N2 

-• 59300E+03 

28.0134 

N 

2 

N2H4 

•22750E+05 

32.0454 

N 

H 

2  4 

AKAR2 

-.25500E+05 

226.2790 

C 

H  N 

0 

14  14 

2 

1 

IPAN 

-.96400E+05 

122.1250 

C 

H  0 

N 

3  10 

3 

2 

HN0  3 

-• 41610E+05 

63.0129 

H 

N  0 

1  1 

3 

PBO 

-*•  10350E+06 

223.1994 

PB  0 

1  1 

DEGDN 

-•  10350E+06 

196.1172 

C 

H  N 

0 

4  8 

2 

7 

AN 

-• 87370E+05 

80*0432 

N 

H  0 

2  4 

3 

MC 

-.17500E+P5 

240.3058 

C 

H  N 

0 

15  16 

2 

1 

HAN 

-• 87600E+05 

96.0430 

N 

H  0 

2  4 

4 

TATB 

-• 33400E+05 

222.1170 

C 

H  N 

0 

3  6 

6 

6 

GLY 

-•15978E+06 

92.0944 

C 

H  0 

3  8 

3 

EGLY 

-• 10873E+06 

62.0688 

C 

H  0 

2  6 

2 

78 


BLAKE 

NAME 

HEAT  OF 

formation 

cal/mole 

MOL. 

WEIGHT 

formula 

CEO  AC 

-.32300E+08 

24265.0000 

C  H  0 

983  1383  692 

02 

-.59300E+03 

31.9988 

0 

2 

CETAC 

-.37340E+08 

28376.0000 

C  H  0 

1179  1579  789 

TEGDN 

-.13850E+06 

240.17IP 

C  H  N  0 

6  12  2  8 

BAO 

-.13100E+06 

153.3394 

BA  0 

1  1 

METRIO 

-.10600E+06 

255.1420 

C  H  N  0 

5  9  3  9 

ETHER 

-.66750E+05 

74.1234 

C  H  0 

4  10  1 

TMAN 

-.74100E+05 

122.1246 

C  H  N  0 

3  10  2  3 

PB2C04 

-.21950E+P6 

490.4086 

PB  C  0 

2  14 

OECA 

-.55120E+05 

138.2522 

C  H 

10  18 

EDDN 

15620E+06 

186.1242 

C  H  N  0 

2  10  4  6 

HN 

-.60130E+05 

95.0578 

N  H  0 

3  5  3 

UDHH 

-• 12690E+Q5 

60.0986 

C 

H  N 

2  8  2 

EOAN 

12360E+06 

124.0962 

C 

H  N  0 

2  8  2  4 

TRIAC 

28018E+06 

218.2070 

c 

H  0 

9  14  6 

HMX 

• 17900E+05 

296.1560 

c 

HON 

4  8  8  8 

NCI  100 

-.18086E+09 

250754.5000 

c 

HON 

6000  8031  8939  1969 

BLAKE 

NAME 

HEAT  OF 
FORMATION 

CAL /MOLE 

MOL* 

WEIGHT 

FORMULA 

NCilOl 

-.18079E+09 

250879.1000 

CHON 

6000  8028  8944  1972 

NC1102 

-.18072E+09 

251003,9000 

CHON 

6000  8025  8950  1975 

NC1103 

-.18066E+09 

251128,8000 

CHON 

6000  8022  8955  1978 

NCI 104 

-• 18059E+09 

251253,8000 

CHON 

6000  8020  8961  1980 

NC1105 

-• 18052E+09 

251378,9000 

CHON 

6000  8017  8966  1983 

NCI 106 

-.18P45E+09 

251504,2000 

CHON 

6000  8014  8972  1986 

NC1107 

-• 18038E+Q9 

251629.6000 

CHON 

6000  8011  8977  1989 

NCllpB 

-• 18p32E+09 

251755,1000 

CHON 

6000  8008  8983  1992 

NC1109 

-• 18025E+09 

251880.8000 

CHON 

6000  8006  8989  1994 

NC111C 

-.I8018E+C9 

252006,5000 

CHON 

6000  8003  8994  1997 

NC1111 

-• 18011E +09 

252132,4000 

CHON 

6000  8000  9000  2000 

NC1112 

-.18005E+09 

252258.4000 

CHON 

6000  7997  9005  2003 

NC1113 

-• 17998E+09 

252384.6000 

CHON 

6000  7995  9011  2005 

NC1U4 

-.17991E+09 

252510.9000 

CHON 

6000  7992  9017  2008 

NC1115 

-.17984E+09 

252637.2000 

CHON 

6000  7989  9022  2011 

NC1116 

-.X7977E409 

252763.8000 

CHON 

6000  7986  9028  2014 

NC1117 

-.17970E+09 

252890.4000 

CHON 

6000  7983  9033  2017 

80 


BLAKE 

NAME 

HEAT  OF 

FORMATION 

CAL/MOLE 

MOL. 

WEIGHT 

FORMULA 

NCU18 

-.17964E+09 

253017.2000 

CHON 

6000  7980  9039  2020 

NC1119 

-.17957E+09 

253144.1000 

CHON 

6000  7978  9045  2022 

NC 1120 

-.17950E+09 

253271.1000 

CHON 

6000  7975  9050  2025 

NC1121 

-.179436+09 

253398.3000 

CHON 

6000  7972  9056  2028 

NC1122 

-,179366+09 

253525.6000 

CHON 

6000  7969  9062  2031 

NC112B 

-.179296+09 

253653.0000 

CHON 

6000  7966  9067  2034 

NC1124 

-.179226+09 

253780.5000 

CHON 

6000  7963  9073  2037 

NC1125 

-.179156+09 

253908. 20p0 

CHON 

6000  7961  9079  2039 

NC1126 

-.179086+09 

254036.0000 

CHON 

6000  7958  9084  2042 

NC1127 

-•  179,026+09 

254163.90,00 

CHON 

6000  7955  9090  2045 

NCI 128 

-.178956+09 

254292.0000 

CHON 

6000  7952  9096  2048 

NC1129 

-.178886+09 

254420.2000 

CHON 

6000  7949  9101  2051 

NC1130 

-.17881E+09 

254548.5000 

CHON 

6000  7946  9107  2054 

NC1131 

-• 17874E+09 

254676.9000 

CHON 

6000  7944  9113  2056 

NC1132 

-.178676+09 

254805.5000 

CHON 

6000  7941  9119  2059 

NC1133 

-.17860E+09 

254934.2000 

CHON 

6000  7938  9124  2062 

NC1134 

-.178536+09 

255063.1000 

CHON 

6000  7935  9130  2065 

81 


BLAKE 

NAME 

HEAT  OF 

FORMATION 

CAL/MOLE 

MOL. 

WEIGHT 

FORMULA 

NCU35 

-• 17846E+09 

255192.0000 

CHON 

6000  7932  9136  2068 

NC1136 

-• 17839E+09 

255321.1000 

CHON 

6000  7929  9142  2071 

NCU37 

-• 17832E+09 

255450.3000 

CHON 

6000  7926  9147  2074 

NCI 138 

-• 17825E+09 

255579.7000 

CHON 

6000  7923  9153  2077 

NCI 139 

-• 17816E+09 

255709.2000 

CHON 

6000  7921  9159  2079 

NC1140 

-•  17811E+09 

255838.8000 

CHON 

6000  7918  9165  2082 

NC1141 

-• 17804E+09 

255968.6000 

CHON 

6000  7915  9170  2085 

NC1142 

-.17797E  +  09 

256098.4000 

CHON 

6000  7912  9176  2088 

NC1143 

-•  17790E+09 

256228.5000 

CHON 

6000  7909  9182  2091 

NC1144 

-.17783E+09 

256358.6000 

CHON 

6000  7906  9188  2094 

NC1145 

— « 17776E+09 

256486.9000 

CHON 

6000  7903  9193  2097 

NCI 146 

-• 17769E  +09 

256619. 300P 

CHON 

6000  7900  9199  2100 

NC1147 

— . 17762E+09 

256749.8000 

CHON 

6000  7897  9205  2103 

NC1148 

-•17755E+09 

25688p.5OO0 

CHON 

6000  7895  9211  2105 

NC1149 

-• 17748E+09 

257011.3000 

CHON 

6000  7892  9217  2108 

NCI 150 

-*•  1Y741E  +09 

257142.3000 

CHON 

6000  7889  9222  2111 

NC1151 

-• 17734E+09 

257273.3000 

CHON 

6000  7886  9228  2114 

82 


BLAKE 

NAME 

HEAT  GF 

FORMATION 

CAL/MOLE 

MOL. 

WEIGHT 

FORMULA 

NCI 152 

-.17726E+09 

257404.6000 

C  4  0  N 

6000  7883  9234  2117 

NC1153 

-• 17719E+09 

257535.9000 

CHON 

6000  7880  9240  2120 

NC1154 

-.17712E+09 

257667.4000 

CHON 

6000  7877  9246  2123 

NC1155 

-.17705E+09 

257799.0000 

CHON 

6000  7874  9252  2126 

NC1156 

-.17698E+G9 

257930. TOCO 

CHON 

6000  7871  9258  2129 

NC1157 

-• 17691E+09 

258062.6000 

CHON 

6000  7868  9263  2132 

NC1158 

-.17684E+09 

258194.6000 

CHON 

6000  7365  9269  2135 

NC1159 

-.17677E+09 

258326.8000 

CHON 

6000  7862  9275  2138 

NC1160 

-.17669E+09 

258459.1000 

CHON 

6000  7860  9281  2140 

NC1161 

-• 17662E+09 

258591.5000 

CHON 

6000  7857  9287  2143 

NC1162 

-• 17655E+09 

258724.1000 

CHON 

6000  7854  9293  2146 

NCI 163 

-» 17648E+09 

258356.8000 

CHON 

6000  7851  9299  2149 

NC1164 

-• 17641E+09 

258989.6000 

CHON 

6000  7848  9305  2152 

NC1165 

-.17634E+09 

259122.6000 

CHON 

6000  7845  9310  2155 

NC1166 

-.176268+09 

259255.7000 

CHON 

6000  7842  9316  2158 

NC1167 

-.17619-+09 

259386.9000 

CHON 

6000  7839  9322  2161 

NC1160 

-« 1761 2E+09 

259522.3000 

CHON 

6000  7836  9328  2164 

83 


SLAKE 

NAME 

HEAT  OF 
FORMATION 

CAL /'MOLE 

MOL. 

WEIGHT 

FORMULA 

NCI 169 

-• 17605E+Q9 

25965/5.8000 

CHON 

6000  7833  9334  2167 

NC1170 

-.17598E+09 

259789.4000 

CHON 

6000  7830  9340  2170 

NC1171 

-•  17590E+09 

259923.2000 

CHON 

6000  7827  9346  2173 

NC1172 

-• 17583E+09 

260057.2000 

CHON 

6i'00  7824  9352  2176 

NC1172 

-« 17576E+G9 

260191.2000 

CHON 

6000  7821  9358  2179 

NC1174 

-• 17569E+09 

260325.4000 

CHON 

6000  7818  9364  2182 

NC1175 

-.17561E+09 

260459.8000 

CHON 

6000  7815  9370  2185 

NCI 176 

-.17554E+09 

260594.3000 

CHON 

6000  7812  °376  2188 

NC1177 

-.17547E+09 

260728.9000 

CHON 

6000  7809  9382  2191 

NC1178 

-.17540E+09 

260863.6000 

CHON 

6000  7806  933  8  219<- 

NC1179 

-.17532E+09 

260998.5000 

C  H  0  N 

6000  7803  9394  2197 

NC1180 

-.17525^+09 

261133.6000 

CHON 

6000  7800  9400  2200 

NCI  l'! 

-• 17518E+C7 

261268.8000 

CHON 

6000  7797  9406  2203 

NC1182 

-.17510009 

261404*1000 

CHON 

6000  7794  9412  2206 

NC1183 

-• 175u3E+09 

261539.5000 

CHON 

6000  779i  9418  2209 

NC1184 

-.17496E+C9 

261675.1000 

CHON 

6000  7788  942*«  2212 

NC1185 

-• 17488E+09 

261810.9000 

CHON 

6000  7765  9430  2215 

BLAKE 

NAME 

HEAT  OF 

FORMATION 

cal/mqle 

MOL. 

WEIGHT 

FORMULA 

NC1186 

-.17481E+C9 

261946.8000 

C  H  0  N 

6000  7782  9436  2218 

NC1167 

-.17474E+09 

262082.8000 

C  H  0  N 

6000  7779  9442  2221 

NC1188 

-• 17466E+09 

262216.9000 

CHON 

6000  7776  9448  2224 

NC1189 

-.17459E+09 

262385.2000 

CHON 

6000  7773  9454  2227 

NC1190 

174522+09 

262491.7000 

ch&n 

6000  7770  9460  2230 

NC 1191 

-.17444E+09 

262628.3000 

CHON 

6000  7767  9466  2233 

NC1192 

-.17437E+09 

262765.0000 

C  H  0  N 

6000  7764  9472  2236 

NC1193 

-.17429E+09 

262901.9000 

C  H  0  N 

6000  7761  9478  2239 

NC1194 

-.17422E+09 

263036.9000 

C  H  0  N 

6000  7758  9485  2242 

NCL195 

-.17415E+09 

263176.1000 

CHON 

6000  7755  9491  2245 

NC1196 

-.174072+09 

263313.4000 

C  H  0  N 

6000  7752  9497  2248 

NC1197 

-.17400E+09 

263450.8000 

CHON 

6000  7749  9503  2251 

NC1198 

-•17392E+09 

263588.4000 

C  H  0  N 

6000  7746  9509  2254 

NC1199 

-• 17385E+09 

263726.1000 

C  H  0  N 

6000  7742  9515  2258 

NCI  200 

-.17377E+09 

263864.0000 

C  H  0  N 

6000  7739  9521  2261 

NC1201 

-.17370E+09 

264002.0000 

CHON 

6000  T  36  9527  2264 

NCI 202 

-•17363E+09 

264140.2000 

C  H  0  N 

6000  7733  9533  2267 

85 


BLAKE 

NAME 

HEAT  OF 

FORMATION 

CAI./MOLE 

MOL  * 

WEIGHT 

FORMULA 

NC1203 

-• 17355E+09 

264278. 5Q00 

C  H  Q  N 

6000  7730  9540  2270 

NC1204 

-.17348E+09 

264417.0000 

CHON 

6000  7727  9546  2273 

NC1205 

-.17340E+09 

264555*6000 

C  H  0  N 

6000  7724  9552  2276 

NCI  206 

-• 17333E+09 

264694.3000 

CHON 

6000  7721  9558  2279 

NC1207 

-.17325E+C9 

264833*2000 

CHON 

6000  7718  9564  2282 

NCI 208 

— *  17318E+09 

264972.2000 

CHON 

6000  7715  9570  2285 

NC1209 

-• 17310E+09 

265111*4000 

CHON 

6000  7712  9577  2288 

NC1210 

-• 17303E+09 

265250.7000 

CHON 

6000  7709  9583  229i 

NC1211 

-• 17295E+09 

265390*2000 

CHON 

6000  7705  9589  2295 

NC1212 

-.17287E+Q9 

265529. 80pC 

CHON 

6000  7702  9595  2298 

NC1213 

-.17280E+C9 

265669.6000 

CHON 

6000  7699  9601  2301 

NCI214 

-• 17272E+09 

265809*5000 

CHON 

60|C0  7696  9608  2304 

NC1215 

-.17265E+09 

265949.6000 

CHON 

6000  7693  9614  2307 

HC1216 

-♦ 17257E+09 

266089*8000 

CHON 

6000  7690  9620  2310 

NC1217 

-*• 17250E+09 

266230*2000 

CHON 

6000  7687  9626  2313 

NCI2X6 

-♦17242E+09 

266370.7000 

chon 

6000  7684  9633  2316 

NC1219 

-• 17234E+09 

266511.3000 

CHON 

6000  7681  9639  2319 

86 


BLAKE 

NAME 

) 

HEAT  OF 

formation 

cal/mole 

MOL* 

WEIGHT 

formula 

NCI  220 

17227E+09 

266652*1000 

c  H  o  N 

6000  7677  9645  2323 

NC1221 

-.17219E+09 

266793*1000 

C  H  0  N 

6000  7674  9651  2326 

NCI  222 

-.17212E+09 

266934.2000 

C  H  0  N 

6000  7671  9658  2329 

MC1223 

-•172Q4E  +09 

267075*4000 

C  H  0  N 

6000  7668  9664  2332 

NCI  224 

-.17196E+09 

267216.9000 

CHON 

6000  7665  9670  2335 

NC1225 

-.17189E+09 

267358*4000 

CHON 

6000  7662  9677  2338 

NCI  226 

-.17181E+09 

267500*1000 

C  H  0  N 

6000  7659  9683  2341 

NC1227 

17173E+09 

267642.0000 

C  H  0  N 

6000  7655  9689  2345 

NC1228 

-• 17166E+09 

267784.0000 

CHON 

600C  7652  9695  2348 

NC1229 

-• 17158E+09 

267926,1000 

C  H  0  N 

6000  7649  9702  2351 

NC1230 

1715CE+09 

2689 68.4000 

C  H  0  N 

6000  7646  9708  2354 

NC1231 

-.17143E+Q9 

268210.9000 

CHON 

6000  7643  9714  2357 

NCI  232 

-*17135E$09 

268353.500® 

CHON 

6000  7640  9721  2360 

NC1233 

-.17127E+09 

268496.2000 

CHON 

6000  7636  9727  2364 

NC1234 

-.17U9E+09 

268639.2000 

C  H  0  N 

6000  7633  9733  2367 

NCI 23 5 

-•17112E+09 

268762.2000 

CHON 

6000  7630  9740  2370 

HC1236 

—.17104E+09 

268925.4000 

chon 

6000  7627  9746  2373 
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BLAKE 

NAME 

HEAT  OF 

FORMATION 

CAL/HOLE 

MOL. 

HEIGHT 

FORMULA 

NC1237 

-.17096E*09 

269068*8000 

CHON 

6000  7624  9753  2376 

NC1238 

-.17088E+09 

269212.3000 

CHON 

6000  7621  9759  2379 

NC1239 

-.17081E+09 

269356.0000 

CHON 

60C0  7617  9765  2383 

NC1240 

-* 17073E+09 

269499.8000 

CHON 

6000  7614  9772  2386 

NC1241 

-.17065E409 

269643.8000 

CHON 

6000  7611  9778  2389 

NC1242 

-•17057 £+09 

269787.9000 

CHON 

6000  7608  9785  2392 

NCX243 

-• 17050E+09 

269932.2000 

CHON 

6000  7605  9791  2395 

NC1244 

-• 17042E+09 

270076.7000 

CHON 

6000  7601  9797  2399 

NC1245 

-• X7034E+09 

270221.3000 

CHON 

6000  7598  9804  2402 

NCI  246 

-.17026E+09 

270366.0000 

CHON 

6000  7595  9810  2405 

NC1247 

-• X7018E+09 

270510.9000 

CHON 

6000  7592  9817  2408 

NCI  248 

-.17010E+09 

270656.0000 

CHON 

6000  7588  9823  2412 

NC1249 

-• 17003E+09 

270801.2900 

CHON 

6000  7585  9830  2415 

NCI  2  50 

-•  16  99  5  E +09 

270946.6000 

CHON 

6000  ?§ 81  9836  2418 

NC125X 

-.16987E+G9 

271092.1000 

CHON 

6000  7579  9842  2421 

NCI  252 

-•  16979E+09 

271237.8000 

CHON 

6000  7576  9649  2424 

NCI  253 

-• 16971E  +09 

271383.7000 

CHON 

6000  7572  9855  2428 
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BLAKE 

NAME 

HEAT  OF 
formation 

CAL/MOLE 

MOL* 

WEIGHT 

FORMULA 

NC1254 

-.16963E+09 

271529.700© 

CHON 

6000  7569  9862  2431 

NC1255 

"• 16955E+09 

271675.8000 

C  H  0  N 

6000  7566  9868  2434 

NCI  256 

-.16947E+09 

271822.1000 

CHON 

6000  7563  9875  2437 

NC1257 

-• 16940E+09 

271968.6000 

CHON 

6000  7559  9881  2441 

NC125B 

-• 16932E+09 

272115.2000 

CHON 

6000  7556  9888  2444 

NC1259 

-, 16924E+C9 

272262.0000 

CHON 

6000  7553  9894  2447 

NC1260 

-.16916E+09 

272409.0000 

CHON 

6000  7549  9901  2451 

NC1261 

-.16908E+09 

272556*1000 

CHON 

6000  7546  9908  2454 

NC1262 

-.1690QE+09 

272703.4000 

CHON 

6000  7543  9914  2457 

NC1263 

-.16 892 £+09 

272850.8000 

CHON 

6000  7540  9921  2460 

NCI  264 

-.16884E+09 

272998.4000 

CHON 

6000  7536  9927  2464 

NC1265 

-.16876E+09 

273146.1000 

CHON 

6000  7533  9934  2467 

NC1266 

-.16868E+09 

273294,0000 

CHON 

6000  7530  9940  2470 

NC1267 

-.16860E+09 

273442.1000 

CHON 

6000  7527  9947  2473 

NCI  268 

-.16852E+09 

273590,3000 

CHON 

6000  7523  9954  2477 

NC1269 

-» 16844E+09 

273738.7000 

CHON 

6000  7520  9960  2480 

NC127Q 

-,168366+09 

273887.2000 

CHON 

6000  7517  9967  2483 
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BLAKE 

NAME 

HEAT  OF 

FORMATION 

CAL/MOLE 

MOL. 

WEIGHT 

FORMULA 

NC1271 

-• 16828E+Q9 

274036. 00] QC 

CHON 

6000  7513  9973  2487 

NCI  2 72 

-• 16820E+09 

274184.8000 

CHON 

6000  7510  9980  2490 

NC1273 

-• 16812E+Q9 

274333.9000 

CHON 

6000  7507  9987  2493 

NCI 2 74 

-• 16804E+09 

274483.1000 

CHON 

6000  7503  9993  2497 

NCI  275 

-.16796E+Q9 

274632.4000 

C  H  0  N 

6000  7500  10000  2500 

NCI 2 76 

-•  16788E  +09 

274781.9000 

c  H  0  N 

6000  7497  10006  2503 

NC1277 

-» 16779E+09 

274931.6000 

C  H  0  N 

6000  7493  10013  2507 

NC1278 

-.16771E+C9 

275081.5000 

C  H  0  N 

6000  7490  10020  2510 

NCI  279 

-• 16763E+09 

275231.5000 

C  H  0  N 

6000  7487  10026  2513 

NC1280 

-.16755E+09 

275381.6000 

C  H  0  N 

6000  7483  10033  2517 

NC1281 

-« 16747E  *-09 

275532.0000 

C  H  0  N 

6000  7480  10040  2520 

NC1282 

-.16739E+09 

275682.5000 

C  H  0  N 

6000  7477  10047  2523 

NC1283 

-• 16731E  +09 

275833.1000 

C  H  0  N 

6000  7473  10053  2527 

NC1284 

-• I6723E  +09 

275984.0000 

C  H  0  N 

6000  7470  10060  2530 

NC1285 

-.16714E+09 

276135.0000 

C  H  0  N 

6000  7467  10067  2533 

NC1286 

-• 16706E+09 

276286.1000 

C  H  Q  N 

6000  7463V 10073  2537 

NC1287 

—  • 16698E  +09 

276437.5000 

C  H  0  N 

6000  7460  10080  2540 
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BLAKE 

NAME 

HEAT  OF 

FORMATION 

CAL/MOLE 

MOL. 

WEIGHT 

FORMULA 

NC1288 

-• 16690E+09 

276588.9000 

C  H  0  N 

6000  7457  10087  2543 

NC1289 

-• 16882E+09 

276740.6000 

C  H  0  N 

6000  7453  10094  2547 

NC1290 

-• 16674E+09 

276892.4000 

C  H  0  N 

6000  7450  10100  2550 

NC1291 

-.16665E+09 

277044.4000 

C  H  0  N 

6000  7446  10107  2554 

NCI  292 

-• 16657E+09 

277196.6000 

C  H  0  N 

6000  7443  10114  2557 

NC1293 

-•16649E+09 

277348.9000 

C  H  0  N 

6000  7440  10121  2560 

NC1294 

-•16641E+09 

277501.4000 

C  H  0  N 

6000  7436  10127  2564 

NC1295 

-.16632E+09 

277654.1000 

C  4  0  N 

6000  7433  10134  2567 

NC1296 

-• 16624E+09 

277806.9000 

C  4  0  N 

60 PP  7430  10141  257C 

NC1297 

-.16616E+09 

277959.9000 

C  4  0  N 

6000  7426  10148  2574 

NC1298 

16608E+09 

2781X3.0000 

C  H  0  N 

6000  7423  10155  2577 

NCi.299 

-.16599E+09 

278266.4000 

C  H  0  N 

6000  7419  10161  2581 

NCI  300 

-• 16591E+09 

278419.9000 

C  4  Q  N 

6000  7416  10168  2584 

NC1301 

-.X6583E+09 

278573.6000 

C  H  0  N 

6000  7412  10175  2588 

NCI  302 

-•  16574E+09 

278727.4000 

C  H  0  N 

6000  7409  10182  2591 

NC1303 

-.X6566E+09 

278881.4000 

C  H  0  N 

6000  7406  10189  2594 

NC1304 

-.16558E+09 

27903&. 6000 

C  H  0  N 

6000  7402  10196  2598 

Si 


8LAKE 

NAME 

HEAT  OF 

FORMATION 

CAL/MOLE 

MOL. 

WEIGHT 

FORMULA 

NC1305 

-.16549E+09 

279190.00 00 

C  4  0  N 

6000  7399  10202  2601 

NC1306 

-.16541E+09 

279344.5000 

C  H  0  N 

6000  7395  10209  2605 

NC1307 

-.16533E+09 

279499.2000 

C  H  0  N 

6000  7392  10216  2608 

NCiaos 

-.16524E+09 

279654.0000 

C  H  G  N 

6000  7388  10223  2612 

N.C1209 

-• 16516E+09 

279809.1*000 

CHON 

6000  7385  10230  2615 

NC13HIC 

-.165C8E+09 

279964*3000 

CHON 

6000  7382  10237  2618 

NC13HL1 

16499 E +09 

280119.7000 

C  H  0  N 

6000  7378  10244  2622 

NC1312 

-.16491E+09 

280275.2000 

C  H  0  N 

6000  7375  10251  2625 

NC1313 

-• 16482E+09 

230431.0000 

C  4  0  N 

6000  7371  10258  2629 

NC1314 

-» 16474E+09 

280586. 9ppp 

C  4  0-  N 

6000  7368  10264  2632 

HC1315 

-• 16465E  +09 

280742.9000 

C  H  0  N 

6000  7364  10271  2636 

NC1316 

-.16457E+09 

280899.20PP 

C  H  0  N 

6000  7361  10278  2639 

NC1317 

-• 16449E+09 

281055.6000 

C  H  0  N 

6000  7357  10285  2643 

NCX318 

-• 16440E+09 

281212.20PP 

C  H  0  N 

6000  7354  10292  2646 

NC1319 

-« 16432E+09 

281369.0600 

C  H  0  N 

6000  7350  10299  2650 

NC1320 

-« 1&423E+09 

281525. 90pp 

C  H  0  N 

6000  7347  10306  2653 

NC1321 

-« 16415E+09 

231683.0000 

C  H  0  N 

6000  7343  10313  2657 

92 


BLAKE 

NAME 

HEAT  OF 
formation 

CAL/MOLE 

MOL, 

WEIGHT 

FORMULA 

NC1322 

-. 16406E  +09 

281840.3000 

chon 

6000  7340  10320  2660 

NC1323 

-.16393E+09 

281997.8000 

C  H  G  N 

6000  7336  10327  2664 

NC1324 

-•16389E+09 

282155.5006 

C  H  0  N 

6000  7333  1«334  2667 

NCI 325 

-.163S1E+09 

282213*3000 

C  H  0  N 

6000  7329  10341  2671 

NCI 3 26 

-.16372E+09 

282471-^3000 

C  H  0  N 

6000  7326  10348  2674* 

NCI 3 27 

16364E+09 

282629  .--S-CffcO 

C  H  0  N 

6000  7322  10355  2678 

NC1328 

-.16355E+09 

282737.8000 

chon 

6000  7319  10362  £681 

NC1329 

16346E+09 

232-9  46.4000 

C  H  0  N 

6000  7315  10369  2 683 

NC1330 

-.16338E+09 

283i05cl000 

CHON 

6S00  7312  1037*  268  8 

NC1331 

-.16329E+09 

283264.0000' 

C  H  0  N 

600)0  7  30  8  10333  2692 

NC1332 

-• 16321E+09 

283423.0000 

C  H  0  N 

6000  7305  10331  2695 

NC1333 

-.16312E+09 

283582.3000 

C  H  0  N 

6000  7301  12398  2699 

NC1334 

-. 16303E+09 

283741,7000 

C  H  0  N 

6000  7298  10405  2<#2 

NC1335 

-.16295E+09 

283901 ,3*000 

C  H  0  N 

6000  7294  10412  2706 

NC1336 

-• 16286E+Q9 

284061.1000 

C  H  0  N 

6000  7291  10419  2709 

NC1337 

-.16278E+09 

284221.1000 

C  H  G  N 

6000  7287  10426  2713 

NC1338 

-.16269E+09 

284381.2000 

C  0  N 

6000  7233  10433  2717 
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BLAKE 

NAME 

MEAT  OF 

FORMATION 

CAL/MOLE 

MOL  . 

WEIGHT 

formula 

NCI  339 

-• 16260E  +09 

284541.5000 

C  H  0  N 

6000  7280  10440  2720 

NC1340 

-» 15252E+09 

284702.0000 

C  H  0  N 

6000  7276  10447  2724 

NC1341 

16243E  +09 

284862*7000 

C  H  0  N 

6000  7273  10455  2727 

NC1342 

-. 16234E  +09 

235023.6000 

C  H  0  N 

6000  7269  10462  2731 

NCI  343 

16225E+09 

285184.7000 

*  H  0  N 

6000  7266  10469  2734 

NC1344 

-• 16217E  +09 

285345.9000 

C  H  0  N 

6000  7262  10476  2738 

NC1345 

-.16208E+09 

235507.30,00 

C  H  0  N 

6000  7259  10483  2742 

NC1346 

-• 16199E+09 

285668.9000 

C  *  0  N 

6000  7255  10490  2745 

NC1347 

-.16 19 It +09 

285830.7000 

C  H  0  N 

6000  7251  10498  2749 

NC 13-48 

-„16182t+09 

285992.6000 

C  H  0  N 

6000  7248  10505  2752 

NC1349 

-  16173E+09 

286154.8000 

C  4  0  N 

6000  7344  10512  2756 

NC1350 

-.16164E+09 

2363  17. 1000 

C  H  0  N 

6000  7240  10519  276C 

NCI  3 51 

16156E  +09 

2864-79.6000 

CHON 

6000  7237  10526  2763 

NC1352 

16147E+09 

2366*42. 3D  00 

C  H  G  N 

6000  7333  10534  2767 

NC1353 

16138E+09 

266BS05.2000 

C  H  0  N 

6000  7230  10541  2770 

NCI 3 54 

16129E+09 

286963.3000 

C  H  0  N 

6000  7226  105  48  2774 

NCI  355 

-.  16120E  +  09 

287131.600© 

C  H  0  N 

6000  7222  10555  2778 

BLAKE 

NAME 

HEAT  OF 

FORMATION 

CAL/MOLE 

MOL. 

WEIGHT 

FORMULA 

NC1356 

-•  16111E +09 

287295.0000 

C  H  0  N 

6000  7219  10563  2781 

NC1357 

-.16103E+C9 

287458. 6000 

C  H  0  N 

6000  7215  10570  2785 

NC1358 

-. 16094E  +09 

287622.5000 

C  H  0  N 

6000  7211  10577  2789 

NC 1359 

-.160856+09 

287786.5000 

C  4  0  N 

6000  7208  10584  2792 

NC1360 

-• 16076E+09 

287950.6000 

C  H  0  N 

6000  7204  10592  2796 

NC1361 

-.160676+09 

288115.0000 

C  H  0  N 

6000  7200  10599  2800 

NC1362 

-.16058E+09 

288279.6000 

CHON 

6000  7197  10606  2803 

NC1363 

-.16049E+09 

288444.3000 

C  H  0  N 

6000  7193  10614  2807 

NC1364 

-.16040E+09 

288609.3000 

C  H  0  N 

6000  7189  10621  2811 

NC1365 

-.160326+09 

288774.4000 

C  H  0  N 

6000  7186  10628  2814 

NCI 3 66 

-.160236+09 

288939.7000 

C  H  0  N 

6000  7182  10636  2818 

NC1367 

•».  16014E+Q9 

289105.2000 

C  H  0  N 

6000  7178  10643  2822 

NC1368 

-• 16005E+09 

289270.9000 

C  H  0  N 

6000  7175  10650  2825 

NCI  3 69 

-• 15996E+09 

289436.8000 

C  H  0  N 

6000  7171  10658  2829 

NCJ.370 

-.15987E+C9 

289602. 90C0 

C  H  0  N 

6000  7167  10665  2833 

NC1371 

-« 15978E+09 

289769,2000 

C  H  0  N 

6000  7164  10673  2836 

NC1372 

-.159696+09 

289935.6000 

C  H  0  N 

6000  7160  10680  2840 

95 


BLAKE 

NAME 

HEAT  OF 

FORMATION 

CAL/MOLE 

MOL, 

WEIGHT 

FORMULA 

NC1373 

-.15960E+09 

290102*3000 

C  H  0  N 

6000  7156  10687  2844 

NC1374 

-.15951E+09 

290269.1000 

C  H  0  N 

6000  7153  10695  2847 

NC1375 

-.15942E+09 

290436*2000 

C  H  0  N 

6000  7149  10702  2851 

NC1376 

-.15933E+09 

290603.4000 

C  H  0  N 

6000  7145  10710  2855 

NC1377 

15924E+09 

290770*8000 

CHON 

6000  7141  10717  2859 

NC1378 

-« 15915E+09 

290938.4000 

C  H  0  N 

6000  7138  10725  2862 

NCI 3 79 

-•  15906E  +  09 

291106.2000 

C  H  0  N 

6000  7134  10732  2866 

NC138P 

-• 15896E+09 

291274*2000 

C  H  0  N 

6000  7130  10740  2870 

NC1381 

-« 15887E  +09 

291442.4000 

C  H  0  N 

6000  7127  10747  2873 

NC1382 

-.15878E+09 

291610.9PPP 

CHON 

6000  7123  10754  2877 

NC1383 

-.15869E+09 

291779.4000 

C  4  0  N 

6000  7119  10762  2881 

NCI 38 A 

15860E+09 

291948. 20PP 

C  H  0  N 

6000  7115  10769  2885 

NC1385 

15851E  +09 

292117. 2000 

CHON 

6000  7112  10777  2888 

NC1386 

-. 15842E  +09 

292286.3000 

C  H  0  N 

6000  7108  10785  2892 

NC1387 

15833E+09 

292455.7000 

C  H  0  N 

6000  7104  10792  2896 

NC1388 

-« 15823E+09 

292625.3000 

C  H  0  N 

6000  7100  10800  2900 

NC1389 

15814E+09 

292795*0000 

C  H  0  N 

6000  7096  10807  2904 
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BLAKE 

NAME 

HEAT  OF 

formation 

CAL/MOLE 

MOL. 

WEIGHT 

FORMULA 

NC1390 

-.158O5E+09 

292965.0000 

C  H  0  N 

6000  7093  10815  2907 

NC1391 

15796E+09 

293135.1000 

C  H  0  N 

6000  7089  10822  2911 

NC1392 

-.15787E+09 

293305.5000 

CHON 

6000  7085  10830  2915 

NC1393 

-« 15778E+09 

293476.0000 

C  H  0  N 

6000  7Q81  10837  2919 

NCI 3V A 

-.15768E+09 

293646.8000 

C  H  0  N 

6000  7078  10845  2922 

NC1395 

-• 15759E+09 

293817.7000 

C  H  0  N 

6CWP  7074  10853  2926 

NC1396 

-.15750E+09 

293988.9000 

C  H  C  N 

6000  7070  10860  2930 

NC1397 

-* 15741E+09 

294160.2000 

C  H  0  N 

6000  7066  10868  2934 

NC1398 

-.15731E+09 

294331.8000 

C  H  0  N 

6000  7062  10875  2938 

NCL399 

15722E+09 

294503.5000 

C  H  0  N 

6000  7058  10883  2942 

NC1400 

-« 15713E+09 

294675.5000 

C  H  0  N 

6000  7055  10891  2945 

appendix  c 


OUTPUT  FROM  TEST  CASE” 
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**  PROGRAM  SLAKE*  VERSION  209.0  ** 


RHO/L  TEMP  PRESSURE  IMPETUS  MOL  WT  CO-VOL  FROZEN  CP(FR)  B (T )  C<TI  GAS  VOL  S  H  E  ADEXP  PHI 

S/CC  K  PSI  FT-L8/L8  GAS  CU  IN  GAMMA  CAL/M-K  CU  IN  IN**6  CU  IN/L8  GIBBS  CAL/HOL  CAL/MOL 

.0373  2189*  1564.  93531.  32.337  24.25  1.2366  10.42  1.869  2.59  730.972  1.43  -816.1  -885.2  1.1325  1.0183 
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No.  Of 

Copies  Organization 

12  Administrator 

Defense  Technical  Info  Center 
ATTN:  DTIC-DDA 
Cameron  Station 
Alexandria,  VA  22314 

1  Office  of  the  Under  Secretary 

of  Defense 

Research  &  Engineering 
ATTN:  R.  Thorkildsen 
Washington,  DC  20301 

1  HQDA/SAUS-OR,  D.  Hardison 

Washington,  DC  20301 

1  HQD A/ DAMA-Z A 
Washington,  DC  20310 

2  HQDA,  DAMA-CSM,  LTC  A.  German 

E.  lippi 

Washington,  DC  20310 

1  HQDA/ SARDA 

Washington,  DC  20310 

1  Commandant 

US  Army  War  College 
ATTN:  Ltbrary-FF229 
Carlisle  Barracks,  PA  17013 

1  Ballistic  Missile  Defense 

Advanced  Technology  Center 
P.  0.  Box  1500 
Huntsville,  AL  35804 

1  Chairman 

DOD  Explosives  Safety  Board 

Room  856-C 

Hof  f man  Bldg .  1 

2461  Eisenhower  Avenue 

Alexandria,  VA  22331 

1  Commander 

US  Army  Materiel  Development 
and  Readiness  Command 
ATTN:  DRCDMD-ST 
5001  Eisenhower  Avenue 
Alexandria,  VA  22333 


No.  Of 

Copies  Organization 

10  Commander 

US  Army  Armament  R&D  Command 
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DRD/.R-LCA 
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L.  Schlosberg 
S.  Westley 

S.  Bernstein 
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C.  Heyman 
Dover,  NJ  07801 

11  Commander 

US  Army  Armament  R&D  Command 
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DRDAR-LCE,  R.  Walker 
N.  Slagg 
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E.  Barrieres 

R.  Davitt 
DRDAR-LCU-CV 

E.  Moore 
C.Mandala 
DRDAR-LCM-E , 

S.  Kaplowitz 
Dover,  NJ  07801 

5  Commander 

US  Army  Armament  R&D  Command 
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